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ABSTRACT 

Weak gravitational lensing is a valuable probe of galaxy formation and cosmology. 
Here we quantify the effects of using photometric redshifts (photo-z) in galaxy-galaxy 
lensing, for both sources and lenses, both for the immediate goal of using galaxies with 
photo-z as lenses in the Sloan Digital Sky Survey (SDSS) and as a demonstration 
of methodology for large, upcoming weak lensing surveys that will by necessity be 
dominated by lens samples with photo-z. We calculate the bias in the lensing mass 
calibration as well as consequences for absolute magnitude (i.e., fc-corrections) and 
stellar mass estimates, for a large sample of SDSS Data Release 8 (DR8) galaxies. 
The redshifts are obtained with the template based photo-z code ZEBRA on the 
SDSS DR8 ugriz photometry. We assemble and characterise the calibration samples 
(~9k spectroscopic redshifts from four surveys) to obtain photometric redshift errors 
and lensing biases corresponding to our full SDSS DR8 lens and source catalogues. 
Our tests of the calibration sample also highlight the impact of observing conditions 
in the imaging survey when the spectroscopic calibration covers a small fraction of 
its footprint; atypical imaging conditions in calibration fields can lead to incorrect 
conclusions regarding the photo-z of the full survey. 

For the SDSS DR8 catalogue, we find crAz/(i+z) ~ 0.096 and 0.113 for the lens 
and source catalogues, with flux limits of r = 21 and r = 21.8, respectively. The 
photo-z bias and scatter is a function of photo-z and template types, which we exploit 
to apply photo-z quality cuts. By using photo-z rather than spectroscopy for lenses, 
dim blue galaxies and L* galaxies up to z ~ 0.4 can be used as lenses, thus expanding 
into unexplored areas of parameter space. We also explore the systematic uncertainty 
in the lensing signal calibration when using source photo-z, and both lens and source 
photo-z; given the size of existing training samples, we can constrain the lensing signal 
calibration (and therefore the normalization of the surface mass density) to within 2 
and 4 per cent, respectively. 

Key words: gravitational lensing: weak — methods: data analysis — galaxies: dis- 
tances and redshifts — galaxies: photometry — cosmology: observations 



1 INTRODUCTION 

The current ACDM cosmological model is dominated by 
the unknown dark compon ents of the universe: dark mat- 
ter and dark energy (e.g.. iKomatsu et al.ll201H ). Gravita- 
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tional lensing, the deflection of light from distant source 
galaxies by intervening masses along the line -of-sight (e.g., 
iBartelmann fc Schneiderl200ll : lRefregierll2003r) , has emerged 
as an enormously powerful astrophysical and cosmologi- 
cal probe. It is not only sensitive to dark energy through 
both cosmol ogical distance meas ures and large-scale struc- 
ture growth IjAlbrecht et al.ll200a ) , but is also sensitive to all 
forms of matter, including dark matter. Measurements of the 
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statistical distortion of galaxy shapes (or weak gravitational 
lensing) due to mass along the line-of-sight have been used in 
numerous studies to constrain the cosmological model, the 
theory of gravity, and the connec tion between galaxies and 



dark matter ( e .g., most recently, iMandelbaum et al 
iFu et al.l l2008l : iReves et all I2OI0I : ISchrabback et alf 



I2OO6I : 
20101 ). 



As a result, many more surveys are planned for the next 
two decades with weak lensing as a major science driver: 
KID^, DEfl HSC0, Pan-STARR^, LSStQ, Euclic0, and 

wfirstQ. 

In order to fully access the information encoded in grav- 
itational lensing, redshift information is essential, as the 
conversion from distortions (gravitational shear) to mass de- 
pends on the lens and source redshifts via the critical surface 
density, expressed as 



E, = 



Ds 



AixGil + ZLYDLDLS 



(1) 



in comoving coordinates (in physical coordinates, the ex- 
pression lacks the factor of (1 -I- 2l)~^). Here, Dl and Ds 
are angular diameter distances to the lens and source, and 
Dls is the angular diameter distance between the lens and 
source. Spectroscopic redshifts provide the best accuracy 
in determining Ec, but obtaining them for large, statistical 
samples of both lenses and sources is prohibitively expensive 
in terms of observing time and instrumentation. As a result, 
upcoming surveys will rely heavily upon less accurate photo- 
metric redshifts (photo-z) derived from multi-band imaging. 
Since we require source samples at higher redshift, and as the 
data pushes into the region of dim galaxies with poor pho- 
tometry, the photo-z may worsen even more. Consequently, 
it is of immediate interest to quantify how limited redshift 
accuracy (due to the use of photometric redshifts) propa- 
gates into lensing results. 

Galaxy-galaxy lensing has been used in the past to 
quantify the connection between lens galaxies or clusters 
and their dark matter (DM) halos, in particular the total 
(average) mass profile around galaxies on > 20kpc scales 
and t he DM halo occupation statistics (JHoekstra et al.l 
I2OO5I : iHevmans e\~^ [2006'; 'Ma ndelbaum et all I2OO6I '): can 
constrain the dark matter power spectrum when used in 
combination with the galaxy 2-point correlation function 
llYoo et all 12004 ICacciato et all l2009l : iBaldauf et al] I2OI0I : 
lOguri fc Takadal2011^ ; and can constrain the theory of grav- 
ity when combined with clustering and redshift-space dis- 
tortions IZhang et alll2007l : iReves et al.|[2010 '). This paper 
addresses the use of photometric redshifts for both sources 
and lenses in galaxy-galaxy lensing, in th e context of the 
Sloan Digital Sky Survey Data Release 8 (jEisenstein et al.l 
I2OIII . and references therein; SDSS DR8 hereafter). Previous 
gala:xy-galaxy lensing studies with SDSS hav e been limited 
to lenses with spectroscopic re dshifts (e.g., ISheldon et al.l 
I2OO4I : IMandelbaum et all Bood ) or photo- 2: with atypically 
high accuracy, such as those for Brightest Cluster Galax;ies 



(e.g., ISheldon et al.ll2009l ). Galaxy-galaxy lensing studies 
with other surveys have typically either involved an unusu- 
ally large numb er of passbands yielding exceptionally good 
photo -z (e.g., iKleinheinrich et al.l 120061 : iLeauthaud et al.l 
l201l[ ). or have had limited area coverage and therefore rela- 
tively poor statisti cs compared with the SDSS studies (e.g., 
IParker et ai1l2007l '). 

Additional work must be done to allow galaxy-galaxy 
lensing to achieve its full potential with large, upcoming 
imaging surveys, and to extend to lens galaxy samples that 
lack spectra in SDSS. Lens galaxy samples that lack spectra 
in SDSS tend to be smaller, dimmer galaxies, or galaxies at 
higher redshifts: it would be interesting to extend galaxy- 
galaxy lensing studies to their DM halo masses and envi- 
ronments, including mass and redshift dependence. While 
studies with other surveys have extended into these regimes, 
they have typically involved deep but very narrow space- 
based data with significant cosmic variance ( Hcvmans ct al.l 
[20061 : ILeauthaud et all I2OI1I ) or wider but still rela tively 
noisy ground-based survey data (|Hoekstra et al.ll2005l ). 

In order to address galaxy-galaxy lensing based on 
SDSS DR8 photo-z, we have calculated a new set of photo- 
metric redshifts for the full flux-limited galaxy sample with 
extinction-corrected model r < 21.8. We applied the pub- 
licly available photo-z code Ziirich Extragalactic B ayesian 
Redshift Analyzer (ZEBRApl. [Feldmann etaLll2006h to the 
SDSS ugriz photometry. In this work, we demonstrate the 
photo-z accuracy by comparing against several spectroscopic 
samples. In the course of this process, we identify several 
concerns related to the observing conditions in the imaging 
survey in regions that overlap the calibration survey, which 
will be generally applicable to any upcoming survey. 

We define the criterion for a "better" photo-z as a 
photo-z that gives minimal scatter in the distribution of 
Zspoc — Zphot, and hence the lowest scatter in the biases of any 
redshift-derived quantities, such as the lensing signal or the 
absolute magnitude. Note that we do not aim for the low- 
est bias in photo-z, according to the philosophy that biases 
can be corrected with a representative calibration sample, 
but rather the lowest uncertainty in the biaqj. This criterion 
results in a low uncertainty in the actual science analysis. 
A brief comparison to other publicly available photo-z cata- 
logues is presented. Improvements in the characterisation of 
photo-z for the photometric galaxies can provide a substan- 
tial boost in statistics for current studies of the large scale 
structure. In addition, working with the photo-z to the sur- 
vey limiting magnitude with SDSS photometry will provide 
an ideal test case for future large surveys that rely upon 
photo-z for all lensing calculations (which are, necessarily, 
dominated by galaxies near the flux limit). 

With the improved and extended SDSS photo-z's, we 
consider the application to galaxy- galaxy lensing. Photo-z 
for weak lensing source galaxies in SDSS was investigated 



http : //www . astro-wise . org/pro jects/KIDS/ 

https : //www . dcirkenergysurvey . org/ 

http : //oir . asiaa. sinica . edu . tw/hsc . php 

http: //pan-Starrs . if a.hawaii . edu/public/ 

http : //www . Isst . org/lsst 

http://sci.esa. int/science-e/www/area/index.cfmTf areaid=102 structure in the photo-z error as a function of redshift, luminosity, 

http://wfirst.gsfc.nasa.gov/ or type. 



* |http : //www . exp-astro . phys . ethz . ch/ZEBRA/ 1 
^ Since the scatter defines the uncertainty, we need a represen- 
tative calibration sample that is large enough to calibrate it. For 
deeper surveys for which such a large representative sample may 
not necessarily be available, the definition of 'ideal photometric 
redshifts' may differ, with a preference for those that show little 
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by iMandelbaum et al.1 120081 . Here, we address some addi- 
tional nuances in the procedures from that work to define 
a fair calibration sample and to use it to estimate the lens- 
ing calibration bias. We then extend that formalism to use 
photo-2 for lenses, thus increasing the possible sample of 
lens galaxies by a factor of ~ 40 over the flux-limited MAIN 
spectroscopic sample (r < 17.77) and colour-selected LRGs 
(luminous red galaxies; r < 19.5). With an eye to using 
photo-z for galaxy-galaxy lensing, we then test how the 
bias and scatter in the photometric redshifts for the SDSS 
photometric galaxy sample affect various derived quanti- 
ties (including the lensing signal calibration, and the esti- 
mated luminosities and stellar masses) by direct comparison 
to their spectroscopic red shift (spec-z) counterparts. While 
iKleinheinrich et al.l (|2005l l identify the need for lens redshift 
information over simply using a lens redshift distribution, 
our work is the first detailed demonstration of how to cal- 
ibrate the effect of the lens photo- 2: errors on other lensing 
observables. 

The outline of the paper is as follows: we describe the 
data and calibration subsets in Sec. [21 and test the imaging 
quality in the regions overlapping the calibration samples 
for consistency with the average survey quality in Sec. [S] 
This latter step is crucial for ensuring that measurements 
using the calibration set are representative of the full SDSS 
DR8 sample of interest. We justify our choice of photo-2: 
method in Sec. 2] The photometric redshift accuracy re- 
sults (bias and scatter) are discussed in Sec. [5] From the 
photometry and the photo-z, we also derive estimated abso- 
lute luminosity and stellar mass, described in Sec. 16.11 and 
Sec. 16.21 respectively. Sec. [6] then describes the resulting bi- 
ases of derived quantities for galaxy-galaxy lensing applica- 
tions. The lensing signal calibration is discussed separately 
in Sec. [T] and conclusions are presented in Sec. El We use a 
flat concordance WMAP7 (|Komatsu et al.ll2011i ) cosmology 
(ilm ~ 0.27, h = 0.702) to calculate luminosity distances 
Dl and angular diameter distances Da from redshifts. 



2 DATA 

Here we describe the data sets used for this investigation: 
the SDSS photometric catalogue and the spectroscopic cal- 
ibration sets. 



2.1 SDSS photometry 

The Sloan Digital Sky Survey (SDSS. lYork et al.ll200ol ) is 
a shallow, wide-field survey that imaged 14 555 deg^ of the 
sky to r ~ 22, and followed up rou ghly two million of the 
detected objec t s spe c troscopically jEisenstein et al.l I2001I : 
[Richards etahl I2OO2I : [Strauss ot al. 200^). T he five-band 
(uqri z) SDSS imaging (Fukugita et al. 19961 : ISmith et al.l 
I2OO2I ) was carried o ut by drift-scanning the sky in pho- 
tometric conditions JHogg et al.1 I2OOII : llvezic et all l2004l : 
iPadmanabhan et al.ll200S ) using a specially-designed wide- 
field camera ( Gunn et al.lll998h . All of the data were pro- 
cessed by completely automated pipelines that detect and 
measure photometric p roperties of objects, and astrometri- 
cally calibrate the data l|Lupton et al.ll200ll : |Pier et al.ll2003l : 
iTucker et aLll2006l 'l. The original goals of SDSS were com- 



pleted with its seventh data release (DR7, lAbazaiian et al.l 
l2009h . 

Accurate galaxy colours are needed in order to compute 
reliable photo-z's. We have chosen our SDSS galaxy sample, 
describe d below, from the mos t recent SDSS-HI Data Re- 
lease 8 (jEisenstein et al.M201ll . DR8). The galaxy colours 
obtained are based o n the model magnitudes MODELMAG 
dStoughton et al.ll2002l ). The 5-band fluxes are estimated us- 
ing a single galaxy model (the better fit of exponential or de 
Vaucouleurs) based on r-band imaging, which is then con- 
volved with the band-specific point-spread function (PSF) 
and allowed to vary only in its amplitude in order to es- 
timate the flux in each band. This procedure leads to a 
consistent definition of the magnitudes across all bands, de- 
spite the different PSFs. This method is superior to PSF- 
matching because it does not require convolutions of the 
data (convolutions lead to correlated noise, making estima- 
tion of the fiux uncertainties challenging). A correction for 
galactic extinc t ion w as imposed using the dust maps from 
ISchlegel et al.l (1 19981) and the extinction-to-reddening ra- 
tios from IStoughton et al.l (|2002h : we only use regions with 
r-band extinction Ar < 0.2. This extinction cut is stan- 
dard for many extragalactic observations, but is also par- 
ticularly important here since dust extinction affects the 
relative magnitudes of the different bands, and the mag- 
nitude corrections become less reliable for higher extinc- 
tions. The pho tometry was calibrated us ing ubercalibra- 
tion procedure (jPadmanabhan et al.ll2008l '). to ensure con- 
sistent calibration across the entire survey (within 1 per 
cent), which is also important for photo-« uniformity. Mi- 
nor corrections were applied to the u and z band to cor- 
rect to AB magnitudes (—0.04 and +0.02, respectively)!^ 
for calculating the photometric redshifts. We note that us- 
ing a more recent estimate on the absolute calibration, 
niAB - msDSS = -0.036,0.012,0.010,0.028,0.040 for the 
ugriz banda^^l. does not significantly modify our results. 

Requirements on the overall data quality for a re- 
gion to be used are that the ubercalibration procedure 
ijPadmanabhan et al.ll2008l ) must classify the data as pho- 
tometric (CALIB_STATUS = 1); the r-band PSF FWHM 
must be smaller than 1.8"; and various Photo fiags must 
indicate no major problems with the object detection and 
PSF estimatioro. In the acceptable regions, a total of 8720 
square degrees, additional cuts were imposed when selecting 
galaxies; see Appendix |X] for details. In regions with mul- 
tiple observations, we first chose the observation with the 
best seeing, and then imposed the magnitude cut. 

The photometric (lens) catalogue is a purely flux- 
limited catalogue of galaxies with photo-z information (this 
work). The full SDSS photometric catalogue has r < 
22, but to avoid the galaxies with very noisy photome- 
try near the flux limit (based on previous flndings, e.g. 
IKleinheinrich et al.ll2005l . that accurate lens redshift infor- 
mation is more important than source redshift information), 
we require r < 21 for our lens sample. Galaxy detection is 
also more stable under different observing conditions with 



http : //cas . sdss . org/dr7/en/help/docs/algorithm . asp' 



http : //howdy . physics . nyu . edu/index . php/Kcorrect 



^^ We require s£ IMAGE.STATUS ^ 4, PSP.STATUS = 0, 
PHOTO.STATUS = 0. 
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this magnitude cut, as discussed in Sec. 13.2.21 resulting in 
a relatively uniform lens galaxy number density across the 
survey. There are 28 036 133 objects in the photometric lens 
catalogue of r < 21 (decreased from 64 750 701 in the full 
r < 22 catalogue). 

The source catalogue, in contrast, goes to a depth 
of r < 21.8 and contains additional information about 
the galaxy shapes (Reyes et al. 2011 in prep.). In order 
to derive galaxy shapes without significant systematic er- 
ror, the source catalogue has further quality cuts in order 
to ensure that the galaxies are suffici ently resolved com- 
pared to the PSF. It is an update of the iMandelbaum et al.l 
((20051) source catalogue, with several technical improve- 
ments and with additional area. As for the old catalogue, 
the galaxy shape measurements are obtained using the RE- 
GLENS pipeline, including PSF corre ction done via re- 
Gaussianization (iHirata fc SeliakI 120031 ) and with cuts de- 
signed to avoid various shear calibration biases. These cuts 
include a requirement that the galaxy be well-resolved in 
both the r and i bands, since the lensing measurements use 
an average of the r and i band shapes. There are 43 378 516 
objects in the source catalogue that satisfy the resolution 
requirement (photo-z quality cuts will reduce this further 
by ~ 10 per cent, as described in Sec. I5.1|l . 



2.2 Calibration data sets 

Quantities derived from photo-z's, such as absolute magni- 
tudes and the lensing signal calibration, will be biased due 
to photo-z erroij^. We estimate these biases (and their un- 
certainty) by measuring them on a representative subset of 
the SDSS catalogue for which spectroscopic redshifts are 
available. This calibration set must satisfy the following cri- 
teria: (1) Target selection is based on apparent magnitude 
only, with no colour selection (the latter often is used in 
spectroscopic surveys to select objects of certain redshift 
ranges and/or types) — unless two samples with complemen- 
tary colour selection can be combined to make an effectively 
flux-limited survey. (2) If targeted off of different photom- 
etry than the imaging survey, then the limiting apparent 
magnitude must be somewhat deeper than the photomet- 
ric survey, because increased flux uncertainty at the lim- 
iting magnitude will randomly scatter objects below and 
above the magnitude threshold. (3) The spectroscopic sam- 
ple should not have any redshift failure modes that have 
a strong preference for particular galaxy types (magnitude, 
colour, or redshift). Note that some of these criteria are not 
absolute, in the sense that the analysis we describe could be 
adapted for calibration data sets that do not satisfy them 
perfectly, but it would significantly complicate the analysis 
and result in greater systematic error. 

There are a limited number of surveys with spectro- 
scopic or spectro-photometric redshift determinations that 
meet these criteria. The rest of this section describes each 
survey that we use for this analysis, noting in particular 



^^ As demonstrated in IMandelbaum et al.l ll200a) . even photo- 
2 that are unbiased on average will cause a bias in the lensing 
signal due to the non- negligible scatter. This is a consequence of 
the nonlinear dependence of the lensing critical surface density 
on the lens and source redshifts. 




Figure 1. Footprints of the SDSS which satisfy our photometric 
quality criteria (8720 deg-^ , shaded) and spectroscopic calibration 
subsamples (labeled). Each calibration subsample covers the full 
photometric depth of SDSS (except PRIMUS-stripe82, see text), 
but only covers ~ 0.2 per cent of the SDSS survey area. 



their survey targeting strategy, our quality cut criteria, and 
the resulting number of calibration galaxies. Figure [T] shows 
the footprints of SDSS and calibration subsamples; Table [T] 
summarises the number of galaxies from each survey that 
have been matched to our SDSS catalogue. Further cuts are 
needed so that the calibration and full samples are represen- 
tative of the same populations (Sec.|3}; there is an additional 
photo- z; quality cut (Sec. l5.l|l . These counts are shown in Ta- 
ble[T]in bold. After the full set of cuts, the parent calibration 
sample has 9 631 sources and 8 5 92 lenses; the former consti- 
tutes a substantial increase over IMandelbaum et al.l l)2008l ). 
which used 2 838 galaxies from only two of the calibration 
samples described below. While we have not addressed the 
spectrosc opic failur e rate on the calibration bias estimates, 
IMandelbaum et al.l (2008 . Sec. 5.5) studied the implications 
of a ~ 5 per cent failure rate and found minimal difference 
on the calibration bias. 



2.2.1 DEEP2 EGS 



The DEEP2 Ga.laxy Redshift Survey ( Davis et al 



Davis et al. 



2003; 



20051 . 



Madgwick et~aLl l2003l : ICoil et al.1 l2004l : 
20071) is a spectroscopic survey usin g the DEep Imag ing 
Muhi-Object Spectrograph (DEIMOs! B^aber et al.ll2003l ') on 
the Keck II Telescope at high spectral resolution {R ~ 5000). 
Of the four DEEP2 fields, we use spec-z's from the Ex- 
tended Groth Strip (EGS) field at RA~14hr, in which the 
targets include galaxies of all colours with Rab < 24.1. For 
R > 21.5, the target selection deviates slightly from a purely 
flux- limited sample. As a result, a small fraction (~ 2 per 
cent) of the SDSS-matched targets have been down- weighted 
based on colour; faint objects with 21.5 < Rab < 24.1 and 
colours suggesting a redshift z < 0.75 receive lower targeting 
priority. However, this down-weighting was found to have 
minimal impact on the derived scatter and biases for the 
overall matched sample. Due to saturation issues, no ob- 
jects brighter than Rab = 17.6 were targeted; these missing 
bright objects constitute a small fraction (< 3 per cent) of 
both our lens and source samples. There is a known redshift 
incompleteness which we discuss later in Sec. 13.4! 

DEEP2 provides redshift quality flags, of which we keep 
flags 3 and 4, which correspond to 95 and >99.5 per cent 
repeatability, respectively. The good-quality (q = 3, 4) red- 
shifts constitute ~ 85 per cent of the SDSS-matched cali- 
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Rcdshift 


Fraction 


Counts 


Survey i 


completeness 


secure 


lens 


source 


EGS 


93 


99 


639 
593 


1060 
968 


zCOSMOS 


85 


>99.5 


775 


1120 


pCOSMOS 


>97 


97 


2506 


3644 


COSMOS (union) 






3281 
2959 


4764 
4293 


VVDS 


55 


97 


643 
564 


1132 
527 


PRIMUS-02hr 


93 


98 


(1436) 


(2588) 


. . . with DEEP2 






1508 
1357 


2908 
1646 


PRIMUS-23hr 


90 


98 


(1612) 


(3023) 


. . . with DEEP2 






1677 
1547 


3307 
1737 


PRIMUS-stripe82 


80 


98 


856 
786 


- 


Total 






8592 


9171 



Table 1. A summary of the properties of our redshift calibration 
samples, including the redshift completeness (in per cent), the 
fraction of secure redshifts corresponding to that completeness 
(in per cent), and the original numbers of lens and source galaxies 
available in each one. The counts in bold are the actual number of 
objects used for calibration, after seeing (Sec. 13. 2.31 1 and photo-2 
quality cuts (Sec. 15.11 1 have been applied. 



bration set, and the fraction are similar for both the lens 
and source sample. Over half of the remaining objects have 
their redshifts confirmed by visual inspection, such that the 
SDSS-matched EGS sample is 99 per cent secure, with 92 
per cent completeness for both the lens and source samples 
(Table [1]). Of the DEEP2 EGS spectroscopic galaxies, 639 
and 1060 objects matched to the SDSS photometric cata- 
logue (r < 21) and source catalogue (r < 21.8 with resolu- 
ti on cuts), respectively. The latter number differs from that 
in lMandelbaum et al.l (|2008r ) because of the use of a new re- 
duction of the SDSS data to create a new source catalogue 
(see section [2TT|) . 



2.2.2 zCOSMOS 



The zCOSMOS-bright survey (JLillv et al.l |2007| . l2009h is 
a magnitude-limited spectroscopic survey on the Cosmo- 
logical Evolution Survey (COSMOS) field (Capak ctld] 
[20071 : IScoville et al.ll2007al lbl: iTaniguchi et al.ll2007l) with the 
Visib le Multi-Object Spectrograph (VIMOS. lLeFevre et all 
120031 ) on the 8m European Southern Observatory Very 
Large Telescope (ESO VLT) at intermediate spectral res- 
olution {R ~ 600). The selection is purely flux-limited at 
15.0 < Iab < 22.5. We choose objects with confidence class 
3 and 4; additionally, we include confidence class 9.5 objects, 
which have redshifts determined from a single emission line 
and which are consistent with the (30-band COSMOS) pho- 
tometric redshifts discussed in section 12.2.31 these consti- 
tute ~1 per cent of the SDSS-matched sample. The photo-z 
agreement is necessary to break the degeneracy between Hq 
and [Oil] 3727, when the doublet appears as a single line due 



to fine broadening. Table [T] lists the number of objects from 
this survey available for photo-« calibration. 

Our match ed zCOSMOS calibr a tion s ample is smaller 
than that from iMandelbaum et al.l (|200g) because ~ 1/3 
of the SDSS imagin g in the COSMOS regio n is classified 
as non- photometric IjMandelbaum et al.ll20ld ). resulting in 
abnormally poor photo-z. Since the non-photometric regions 
are not representative of our source catalogue (for which 
we impose a photometricity cut), we eliminate the region of 
COSMOS for which the SDSS data are not photometric from 
consideration in this analysis, leaving 775 and 1120 lens and 
source galaxies (before further cuts that will be described 
below). An additional consideration regarding the COSMOS 
calibration sample is that the sky noise level in the SDSS 
imaging in the COSMOS region is significantly higher than 
what is typical for the SDSS survey overall. As we will show 
in Sec. 13.21 the consequence is an atypical deficit of fainter 
galaxies in the COSMOS calibration sample; we discuss how 
this deficit is handled in Sec. 17.21 



2.2.3 COSMOS photometric redshifts (pCOSMOS) 

The first new calibration sample used in this work is the sam- 
ple of flux-limited non-zCOSMOS galaxies in the CO SMOS 
region, with photometric redshifts (Ilbert et al. 2009, pCOS- 
MOS hereafter). Given the fiux limit of the SDSS catalogue 
with respect to the COSMOS observations, and the accuracy 
demanded in the applications described in this paper, these 
COSMOS photo-z are effectively the same as spectroscopic 
redshift 0. 

The pCOSMOS photo-z are obtained from a x^ 
template-fitting method, Le Pharg^^l. Highly accurate 
photo- 2: result from the use of deep photometry in 30 bands 
(primary bands are u* , Bj, Vj, g'^ , r"*", i^ , z"*", and K), ob- 
served at various telescopes, primarily from Subaru and the 
Canada- France-Hawaii Telescope (CFHT). Comparing the 
COSMOS photo-z with the matching zCOSMOS redshifts, 
we estimate the scatter in the pCOSMOS calibration subset 



to be (Ta 



z/{l + z) 



0.009, with an outlier rate of ~ 1.2 per 



cent for both the lens and the source samples (we find the 
bias to be negligible). COSMOS galaxies that have a zCOS- 
MOS spectroscopic redshift have been removed, so that the 
zCOSMOS and pCOSMOS catalogues are disjoint (though 
they trace similar large-scale structures). There are 2506 and 
3644 pCOSMOS matches to the photometric lens and source 
catalogues, respectively. 



2.2.4 VVDS 

The ne xt new calibration sa mple is the VIMOS VLT Deep 
Survey (|Le Fevre et al.ll2005l . VVDS), another spectroscopic 
survey using the VIMOS instrument on the ESO- VLT. We 
use the public spectroscopic catalogue in the VVDS-F22 
field, part of the VVDS-wide survey. This survey uses a 
lower resolution {R ~ 230 versus 600) but similar exposure 



^* That is, we have verified that increasing the true redshift un- 
certainty by up to a few per cent will not degrade our ability to 
estimate the bias in derived redshift quantities, because the errors 
in the SDSS photo-2 will always dominate. 
^^ http : //www . oamp . f r/people/arnouts/LE_PHARE.html| 
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time (3000 versus 3600 sec) compared to the zCOSMOS- 
bright survey. The VVDS-wide target selection is purely 
magnitude-limited (17.5 < Iab < 22.5) with no colour se- 
lection. We select those galaxies with spectroscopic quality 
flags 3 and 4, with ~96 and ~99 per cent reliability, respec- 
tively. When combined, these constitute only ~55 per cent of 
the SDSS-matched sample, suggesting a relatively low red- 
shift success rate for this sample. Additionally, galaxies with 
qualities 13, 14 (broad emission line) or 23, 24 (serendipi- 
tous observations) are included; however these constitute 
a small fraction (3 per cent) of the SDSS-matched sam- 
ple. In total, we use 643 lens and 1132 source galaxies with 
VVDS spectroscopic redshifts for SDSS photo-z calibration, 
see Table [T] Because of concerns about the high spectro- 
scopic redshift failure rate, we have (1) checked the colour 
and magnitude distributions of the VVDS-matched sample 
(Section l3.ip and (2) done all of the calculations of biases due 
to photo-2 error separately for individual spectroscopic sub- 
samples. We find (Section VJ.2\ no systematic tendency for 
the VVDS results to differ with any of the others, suggesting 
that the redshift failure rate has not excessively biased the 
calibration sample properties. 

2.2.5 PRIMUS 

The final additions to our calibration sa mple come from 
the PRIsm MUlti-object Survey (PRIMUS : ICoil et al.|[201ol . 
Cool et al. 2011 in prep.), a newly-completed survey that ob- 
tained spectro-photometric redshifts for ~ 120 000 galaxies 
to i = 23 over an area of 9.1 deg^. On the Inamori Magellan 
Areal Camera and Spectrograph (IMACS) on the Magel- 
lan I Telescope, PRIMUS uses a low-dispersion prism with 
7? ~ 40 to efficiently survey wide areas and achieve a red- 
shift accuracy of o"z/(i+z) ^ 0.005 out to z = 1.2. From the 
PRIMUS low-resolution spectra, the general shape of the 
spectrum, in addition to emission or absorption lines, is used 
to infer the redshift (Cool et al. 2011 in prep.). PRIMUS is 
generally a flux-limited survey (no colour selection); how- 
ever there are a few relevant exceptions to this rule that 
we will discuss shortly. The target selection in most tar- 
get fields used full and sparse sampling for i < 22.5 and 
22.5 < i < 23.5, respectively, although the limiting magni- 
tudes and reference band for these categories vary depending 
on the target field. 

We included data from 3 fields, which we will denote 
PRIMUS-02hr, PRIMUS-23hr, and PRIMUS-stripe82. The 
PRIMUS-02hr and 23hr fields overlap the DEEP2 02hr and 
23hr spectroscopic survey fields (these DEEP2 fields are 
BRI colour-selected to include almost exclusively z > 0.7 
galaxies). The magnitude limits in these fields for full and 
sparse sampling were R < 22.8 and 22.8 < R < 23.5, re- 
spectively. The primary targeting strategy in these PRIMUS 
fields was to target the complement of the DEEP2 spectro- 
scopic survey; hence the union of the PRIMUS and DEEP2 
surveys in the 02hr and 23hr fields is essentially a flux- 
limited sample (|Coil et al.ll2010l ). We have verifled that the 
relative weighting of DEEP2 to PRIMUS, which accounts 
for the targeting efficiency and survey area differences, is 
very near unity in the SDSS-matched sample in both fields. 
A fraction of the PRIMUS targets also have high-quality 
redshifts in DEEP2; for these, the spectroscopic redshifts 
from DEEP2 were used. The supplements from the DEEP2 



survey with quality flag >= 3 constitute ~5 and ~10 per 
cent of photometric lens (r < 21) and source (r < 21.8) 
catalogues, respectively. We limit our PRIMUS redshifts to 
those having quality flag of 3 or 4. The corresponding relia- 
bilities are estimated at 93 and 97 per cent, respectively. 

Like the other two flelds, the PRIMUS-stripe82 cali- 
bration field overlaps the SDSS stripe-82 region, but in this 
case the PRIMUS targeting was carried out from single- 
epoch SDSS imaging to r < 22 or i < 21 with no colour 
selection. The fact that this targeting was carried out using 
a different (earlier) reduction of the SDSS imaging, in one 
particular observing run, results in targeting incompleteness 
at the faint end, with respect to the DR8 photometry. This 
issue, which is caused by the scatter in the derived magni- 
tudes between different reductions, is apparent in Figure [5] 
below. Due to this incompleteness, which results in this re- 
gion being quite different from the others with respect to the 
abundance of r > 21 galaxies, we have opted to drop the ~lk 
objects for the dimmer (source) catalogue. The matches to 
the brighter lens catalogue were used for our calibration; 
these include 856 objects (80 per cent completeness). The 
combination of the three PRIMUS fields adds over 3000 lens 
and 5000 source calibration galaxies to our calibration sam- 
ple. 

Table [T] shows a summary of the number of galax- 
ies passing the cuts from each of the above calibration 
fields, and the remaining fraction after cuts described be- 
low (Sec. 13.2.31 Sec. 15. 1[) are applied. 



3 CALIBRATION SET ADJUSTMENTS 

The spectroscopic calibration sample is an extremely small 
subset of the SDSS photometric sample (Fig. [l]). To ensure 
that conclusions derived from these calibration samples are 
applicable to the full SDSS sample, the following differences 
between calibration subsamples must be addressed: 

• The differing spectroscopic survey completeness in each 
subsample. 

• The differing SDSS observational conditions in the 
imaging data overlapping each spectroscopic calibration sub- 
sample (SDSS completeness). 

• The specific large-scale structure (LSS) in each of the 
calibration fields (sample variance). 

We generate two calibration sets, one each for the lens and 
source catalogue. The first two issues are most relevant for 
the source catalogue, as its faint galaxies are less robust to 
these issues. 



3.1 Spectroscopic completeness 

The SDSS-matched spectroscopic calibration samples de- 
scribed in Sec. 12.21 have varying degrees of magnitude and 
colour completeness. Here we address the targeting incom- 
pleteness of the spectroscopic set with respect to the SDSS 
sample. The completeness of the SDSS imaging data itself 
is addressed in Sec. 13.21 

Figure[2]shows the r-band magnitude distribution of the 
spectroscopic sample, matched to the source (r < 21.8) and 
photometric (r < 21) SDSS catalogues (the thick and thin 
solid lines, respectively). The corresponding dashed lines 
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Figure 2. The r-band magnitude distribution of the spectro- 
scopic caUbration samples, matched to the source (r < 21.8, thick 
solid line) and photometric (r < 21, thin solid line) SDSS cata- 
logues, respectively. The corresponding dashed lines indicate the 
underlying SDSS distribution in the same area. There is a notice- 
able discrepancy within ~ 0.5 mag of the limiting magnitude for 
the PRIMUS-stripe82 spectroscopic set relative to SDSS, indicat- 
ing spectroscopic incompleteness for the source sample at r > 21 
(see text). 



epoch imaging. The earlier reduction resulted in magnitude 
incompleteness at the faint end, caused by the scatter in 
the derived magnitudes between different reductions. Since 
the PRIMUS-stripe82 source-matched catalogue also shows 
a lack of higfi redshift {z > 0.7) objects which are important 
for source catalogue calibration, we remove this field from 
the source calibration set (but not from the lens calibration 
set, since there are fewer discrepancies there). 

Similarly, the thick and thin solid lines in Figure |3] show 
the g — r colour distribution of the spectroscopic sample. 
Differences in the dashed (underlying SDSS distribution) 
and solid (distribution of the high-quality spectroscopic red- 
shifts) lines indicate colour dependence of the spectroscopic 
redshift success rate. Most samples show no noticeable dis- 
crepancy in the g — r colour distributions. We have examined 
all 4 colours {u — g, g — r, r — i, i — z), and found VVDS shows 
some discrepancy in the g — r colour (shown) , and to a lesser 
degree, in r — i. Although this is alarming, we have other- 
wise not found anomalous behaviour in the VVDS subsam- 
ple; the redshift distribution, the photo- 2; spectral template 
type distribution, and the various derived biases are consis- 
tent with those for the other subsamples. Hence we keep the 
VVDS calibration field. We note that the PRIMUS-stripe82 
source-matched sample shows discrepancy in the r — i and 
i — z colours, but in such a way that mimics the colour dis- 
tribution of the brighter sample; this indicates consistency 
with the incompleteness in the r magnitude distribution. 




).5 0.0 0.5 1.0 1.5 2.0 2.5 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 
g-r coior g-r coior 



Figure 3. The g — r colour distribution of the spectroscopic 
calibration samples, matched to the source (r < 21.8, thick solid 
line) and photometric (r < 21, thin solid line) SDSS catalogues, 
respectively. The corresponding dotted lines are the underlying 
SDSS distribution in the same area. Difference in the dotted and 
solid lines indicates colour dependence of the success rate of the 
spectroscopic redshifts. VVDS shows some colour dependence (g— 
r) of the redshift success rate (see text). 



show the underlying SDSS distribution in the same area. 
The PRIMUS-stripe82 field shows a noticeable discrepancy 
in the number counts at r > 21 between the spectroscopic 
redshift sample and the underlying magnitude distribution. 
The magnitudes used for target selection in this field came 
from an earlier photometric reduction of the SDSS single- 



3.2 SDSS completeness 

Here we address completeness of the SDSS itself. Different 
SDSS fields exhibit different r magnitude distributions, due 
to the varying observational conditions. This can be seen 
in the different panels of Fig. [2] which have different dashed 
histograms (which are for the source samples in these regions 
before requiring a match in the spectroscopic data), though 
in the absence of other information we cannot rule out that 
these variations are due to sampling variance rather than 
observing conditions. Here we will demonstrate how the non- 
uniform seeing and sky noise over the photometric survey 
area, which affects the depth to which an object can be 
detected, can give rise to the differences in this figure. For 
the calibration set to fairly represent the SDSS as a whole, 
we need to understand how the observing conditions in each 
field differ from the median of the whole survey, and correct 
for any severe deviations from typical conditions. 



3.2.1 SDSS observing conditions 

Figure |4] shows the SDSS observation conditions in r-band, 
for the source galaxies in each of the calibration fields and 
in the full SDSS survey. Here we have used galaxies that are 
well-resolved to trace the observing conditions, which means 
that we will be skewed towards better observing conditions 
(as compared to if we had done the calculation using pure 
areas). This methodology explains why the seeing is notice- 
ably better than the typical value that is commonly used for 
SDSS, 1.4". Points indicate a median value, and error bars 
indicate the 68th percentile. To calculate sky noise (horizon- 
tal axis), we approximate the Poisson noise due to the sky 
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Figure 4. SDSS r-band observation conditions (sky noise, as 
described in the text, and seeing) for each of the calibration fields 
and the SDSS full survey. 



and dark current as 



sky 
gain 



(Taky (counts) ~ \\ — ; h dark current 



(2) 



This is a reasonable approximation for galaxies with r > 20, 
where the Poisson noise due to the galaxy flux itself is neg- 
ligible. We then convert this to the sk y noise in nanomag- 
gies u sing the calibration from ubercal ( Padmanabhan et al.l 
I2OO81 ): this noise in nanomaggies determines the S/N for a 
galaxy with a given size and flux, at fixed seeing and galac- 
tic extinction A,-. Then, we multiply the sky noise estimates 



by 10" 



in reality, of course, extinction modulates the 



flux, but this has the same effect on the S/N as increasing 
the sky noise by this factor. Finally, we have divided out by 
the median value of lO'''^ ""fsky for the survey. When the 
sky noise is large, more galaxies at a given flux will fail the 
S/N > 5 object detection filter; this is also true for worse 
seeing, since that filter is imposed within a PSF rather than 
using the entire galaxy flux. 

The r-band PSF FWHM in arcsec is shown on the verti- 
cal axis. In poor seeing, the observed magnitudes are noisier 
(because the galaxy is spread out over more pixels and has 
greater sky noise contribution), and star-galaxy separation 
is more challenging (more galaxies get classified as stars). 

As shown, the observing conditions in the calibration 
fields are not typical of the full SDSS DR8 sample. The 
EGS is within the la range of observing conditions for both 
seeing and sky noise. The COSMOS data have relatively 
high sky noise for more than half the area, which will lead 
to reductions in galaxy S/N of ~ 20 per cent compared 
to that for similar galaxies in the EGS. PRIMUS-02hr and 
PRIMUS-23hr are marginally within the typical region for 
sky noise, and VVDS and PRIMUS-stripe82 are at higher 
sky noise; but all four fields have significantly better see- 
ing than that of the full SDSS. The reason for this trend is 
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Figure 5. An illustration of how observing conditions can modu- 
late the observed galaxy number densities, based on comparisons 
of pairs of runs with identical area coverage but different observ- 
ing conditions in stripe 82. Each row utilizes a different pair of 
runs with conditions given on the right. Black lines are for the 
source catalogue, red lines are for the full photometric catalogue 
(no lensing selection). Left column: Histogram of r-band model 
magnitudes, with the different colours indicating which catalogue 
was used, and different line types indicating the observing condi- 
tions as labeled on the right-hand side. Right column: Solid lines 
show the ratio of the histograms with different observing condi- 
tions from the left column, always putting the dashed histogram 
in the numerator. The dotted lines show the ratio of the total 
number of galaxies in the region as defined using the full sample 
(averaged over magnitudes). 



that all of these samples are on stripe 82, which has many 
observations. Since we select towards observations with bet- 
ter seeing (preferable for lensing) when creating the source 
catalogue, we get atypically good seeing. The fact that the 
extinction is on the higher end for these regions increases 
the effective sky noise, though it is not the only (or even the 
main) reason why these calibration samples have higher sky 
noise than is typical. As we will see in the next section, to 
lowest order we expect that the changes in sample properties 
due to having atypically good seeing will oppose (and even 
outweigh) the changes in sample properties due to atypically 
high sky noise for these stripe 82 samples. 



3.2.2 Effect on number counts 

In Fig. O we illustrate the effect of observing conditions 
on the overall counts and magnitude distribution of galax- 
ies, which is helpful in understanding the behaviour of the 
galaxy counts in the different calibration samples. We make 
this comparison by selecting pairs of runs overlapping the 
PRIMUS-02hr and PRIMUS-23hr regions, with different ob- 
serving conditions but exactly the same area coverage. Thus, 
any differences in the galaxy apparent magnitude distribu- 
tions in the pairs of runs are due to the observing conditions. 



© 0000 RAS, MNRAS 000, 000-000 



Lensing with photometric redshifts 9 



rather than intrinsic differences in the galaxy populations. 
The four rows correspond to four different scenarios, with r 
model magnitudes on the horizontal axis and number counts 
on the vertical axis. At left are the counts (for both fields, 
source and photometric lens catalogues). At right the ra- 
tio between counts for the runs with two different observing 
conditions is shown. In the top row, results are given for 
nearly identical (within 10 per cent) sky noise, but 45 per 
cent difference in the seeing. (This is done by comparing 
runs that are ~ 20 per cent better and worse than the sur- 
vey median seeing.) The second row corresponds to seeing 
that is nearly identical (within 10 per cent) and a bit bet- 
ter than the survey median, but sky noise that differs by 
30 per cent. The third row shows nearly identical conditions 
for the seeing and sky noise, to illustrate the magnitude of 
the differences that can occur just due to Poisson noise. The 
fourth row shows the case of two effects going in the same 
direction: 25 per cent difference in seeing and 15 per cent 
difference in sky noise, where one run is better in both cases 
(better seeing and lower sky noise). 

For the photometric catalogue, the effect of seeing on 
the galaxy number counts is small (< 5 per cent, as shown 
in the top row) and nearly constant over all magnitudes. 
This change in counts may be due predominantly to the 
star-galax;y selection. For the source catalogue, which has 
apparent size selection, the effect of seeing is larger (as we 
would expect), and varies more strongly with magnitude, 
from a nearly constant 20 per cent below r < 21 to typically 
50 per cent above r > 21. 

The signature of sky noise, shown in the second row, 
is somewhat different. It preferentially reduces the num- 
ber counts at fainter magnitudes for the purely flux-limited 
galaxy sample. The overall effect of 5 per cent is actually 
dominated by the fainter parts of the catalogue, with al- 
most no effect below r < 21.4 and a gradually increasing 
effect at fainter magnitudes, up to 35 per cent. In contrast, 
for the source catalogue, the average difference is more sig- 
nificant (15 per cent) and becomes noticeable at brighter 
magnitudes, from r >~ 20.8. The third row shows that for 
identical observing conditions, the samples that are selected 
in both the photometric and source catalogue are statisti- 
cally identical, with < 5 per cent fluctuations in the magni- 
tude histograms. Examination of the actual galaxy samples 
indicates significant noise, with fainter galaxies being scat- 
tered in and out of the sample, such that the fraction of 
galaxies in both catalogues is ~ 70 per cent. Finally, the 
fourth row is not at all surprising considering the first and 
second; we see both the effects of seeing and of sky noise. 

From this analysis we see that the fiux-limited photo- 
metric lens sample, with r < 21, is not strongly affected by 
the changes in observing conditions, because these effects 
are only very significant (> 5 per cent) in objects fainter 
than r ~ 21.4. In contrast, the source number density is 
strongly affected by observing conditions. We now turn to 
the corrections we make to account for these effects in our 
calibration samples, to make them more representative of 
the fuU SDSS area. 



3.2.3 Corrections for observation conditions 

Figure |6] explicitly compares the magnitude distribution for 
each of the source calibration samples to that of the whole 




Figure 6. The r-band magnitude distribution of the source cata- 
logues for each spectroscopic subsample (line histograms), in con- 
trast to the SDSS distribution (grey block histograms). The dif- 
fering SDSS observing conditions cause variable magnitude com- 
pleteness near the limiting magnitude r ~ 21.8. We adjust the 
galaxy population in VVDS, PRIMUS-02hr, and PRIMUS-23hr 
fields to match the SDSS completeness by systematically trim- 
ming the lowest resolution objects (dotted lines), such that the 
resulting r magnitude distribution matches that of SDSS (see 
text). 



of SDSS. Objects in the fields that overlap SDSS stripe 82 
(VVDS, PRIMUS-02hr and PRIMUS-23hr) tend to have 
high completeness near the limiting magnitude; this is be- 
cause there are many repeat observations in this stripe, and 
the source catalogue reductions choose the observation of 
each galaxy that has the best seeing. 

Since stripe 82 only covers a small fraction of the SDSS 
footprint (~3 per cent) and is not representative of the whole 
of SDSS, we would like to correct for the excess faint objects. 
These objects are removed as follows: (1) we bin the galaxies 
by their r magnitudes, and (2) for each magnitude bin, we 
remove the lowest apparent resolution galaxies until the re- 
sulting magnitude histogram matches that from all of SDSS. 
This process removes ~ 40 per cent of the objects in a given 
field. (Figs. [5] and [3] do not include this cut.) 

The COSMOS field has poor sky noise, and in general 
has unusually bad observational conditions compared to the 
rest of SDSS (e.g., a third of the COSMOS field is not clas- 
sified as photometric in SDSS, so we have eliminated that 
portion of the field). This leads to a deficit of faint galaxy de- 
tections; however, we cannot artificially add objects that are 
missing from the catalogues. Hence, when using the COS- 
MOS source calibration set, we have applied empirical cor- 
rections to the derived values (see Sec. 17.21) . 



3.3 Sample variance 

Each individual field has distinct large-scale structure due 
to galaxy clustering. Thus, the observed dA''/d2: in each field 
differs from the true underlying redshift distribution dN/dz, 
which is unknown but is expected to be smooth. We model 
this smooth distribution by fitting to a functional form (our 
bias analysis does not depend significantly on the detailed 
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Figure 7. The observed redshift distributions of the lens caU- 
bration subsamples are shown as thick line histograms. The fits 
to the underlying redshift distribution (thin lines), along with the 
best-fitting parameters, are also shown. 
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Figure 8. Analogous to Fig.[7l except these are for the source cal- 
ibration subsamples, along with the combination of all the source 
calibration subsamples excluding the COSMOS field (see text). 



form of the curve, as discussed in iMandelbaum et al.ll2008n 
AN 
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which has a mean redshift of 



(2) = 



V2 2H.r[(«-H)/2] 
r[a/2] 
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Figures [7] and |8] show the true redshift distributions for 
each of the lens and source calibration samples, along 
wi th their fitted curve s and parameters. Following Sec. 4.2 
of IMandelba um et al.l (J2008), the histograms are x^"fit to 
Eq. ([3} with flat weights while constraining the area under 
the curve to equal to the galaxy count A'gai; the redshift 
histograms of binning width Az = 0.05 were bootstrap re- 
sampled to obtain LSS error estimates to the fitting curve. 



Table 2. Fitting parameters for redshift distributions of the 
subsamples. The source samples for WDS, PRIMUS 02hr and 
PRIMUS 23hr have the seeing cuts (Sec. 13. 2. 31 1 applied; the 
redshift distribution parameters before and after photo-^ cuts 
(Sec. 15.11 1 are listed. 



Sample 



A', 



gal 



(^> 



Photometric (lens) sample with r<21 



EGS 


639 


0.209 


3.050 


0.336 


WDS 


643 


0.224 


2.710 


0.338 


COSMOS 


3281 


0.207 


3.139 


0.339 


PRIMUS 02hr 


1508 


0.243 


2.325 


0.334 


PRIMUS 23hr 


1677 


0.220 


2.987 


0.349 


PRIMUS stripe82 


856 


0.223 


2.718 


0.336 


Photometric (lens) sample 


with r<21 and photo-2 quality cuts 


EGS 


593 


0.204 


3.277 


0.343 


WDS 


564 


0.219 


3.009 


0.350 


COSMOS 


2959 


0.183 


4.106 


0.350 


PRIMUS 02hr 


1357 


0.230 


2.665 


0.343 


PRIMUS 23hr 


1547 


0.209 


3.328 


0.354 


PRIMUS stripe82 


786 


0.213 


3.061 


0.345 


Source sample with 


r<21.8 


(includes seeing 


cuts) 


EGS 


1060 


0.289 


2.399 


0.404 


WDS 


609 


0.253 


2.101 


0.391 


COSMOS 


4764 


0.253 


2.560 


0.368 


PRIMUS 02hr 


1834 


0.344 


1.854 


0.411 


PRIMUS 23hr 


1896 


0.308 


2.257 


0.416 


Source sample with 


r<21.8 and photo-z quality cuts 


EGS 


968 


0.275 


2.689 


0.411 


WDS 


527 


0.293 


2.352 


0.405 


COSMOS 


4293 


0.234 


3.085 


0.380 


PRIMUS 02hr 


1646 


0.329 


2.071 


0.422 


PRIMUS 23hr 


1737 


0.298 


2.454 


0.422 



The distributions shown here were derived after imposing 
cuts to correct for different seeing in the calibration samples 
on stripe 82, as discussed in Sec. 13.2.31 When further photo- 
z quality cuts are applied (Sec. I5.1|l , the distribution shifts 
somewhat; the fit parameter values for both cases are listed 
for each catalog and subsample in Table [21 

The errors in the fits to the redshift distribution for each 
sample come primarily from sampling variance; thus, when 
fitting to the overall calibration sample redshift distribu- 
tion, these fiuctuations are reduced (more significantly than 
one would expect fiuctuations due to Poisson noise to be re- 
duced). The significant smoothing out of LSS can be seen in 
the bottom right panel of Fig. |S] where we plot the sum of 
all source subsamples, excluding the COSMOS region (this 
region was excluded due to its difference in sky noise which 
should result in a different intrinsic redshift distribution). 
This curve is our best estimate for the source catalogue red- 
shift distribution; a similar smoothing in LSS is seen for the 
lens sample, for the sum of all 6 subsamples. These figures 
show that the upper limit in redshifts is around z ~ 0.8 
(z ~ 1.0) for the lens (source) catalogues, with a mean red- 
shift of 0.34 (0.41). We note that the parameters z, and a 
are degenerate to some extent (small a brings the peak to 
higher redshifts, while high z* spreads out the overall distri- 
bution), and hence there is considerably more uncertainty 
in the parameters in [z«,a] than in the mean redshift (z). 
The fit parameter errors can be estimated by modulating 
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the r edshift bin counts to simulate LSS (|Mandelbauni et alj 
120081 ). or by looking at the variation in the fit values of the 
different subsamples. 



3.4 Other incompleteness 

The EGS matched samples have a deficiency of 2 ~ 0.3 
galaxies, which is due to instrumental limitations of the sur- 
veM^^I (Newman et al. 2011 in prep., priv. comm.). Such a de- 
ficiency cannot easily be remedied; however, we find that, at 
this redshift, this deficiency does not affect the derived lens- 
ing signal calibrations when using source photo-z, which is 
more sensitive to the distribution of high-redshift (z > 0.6) 
objects. However, this does affect the calibration for lens 
photo-2, and care must be taken when the EGS lens cali- 
bration subset is used on its own (Sec. 17.3)) . 



4 PHOTO-Z METHOD 



We c hose the template-based ZEBRA code (jFeldmann et al.l 
I2OO61 ) to estimate our photo-z's based upon the SDSS ugriz 
bands. There are several public SDSS photo-2 catalogues 
for the full photometric sam ple, all of wh i ch are based on 
training-set methods: two b ylOvaizu et al.l (120081). available 
in the DR8 CAS data bas(E!|, one bv iBudavari et aL JpOOd ) 
in SDSS D R7 (lAbazaiiaiTet al.lJ2009l ') CAS databasfT 



and 



another bv lCunha et all l|2009l fT" 

Thes e photo-z methods were included among those 
tested in iMandelbaum et al.l l|2008r ). for source galaxy red- 
shifts. However, as those estimates were based on earlier 
Data Releases (5 and 6), including more preliminary photo- 
z estimation procedures and photometr y that lacked uber- 
calibration, the results of the tests in IMandelbaum et al.l 
l|2008r ) do not necessary apply to these new DR7 and DR8 
versions. In Appendix|B] we compare the photo-z's in those 
publicly available SDSS DR8 catalogue, and give our rea- 
sons for choosing to use ZEBRA, another template-based 
method that will be described below. 

Template-fitting photo-z algorithms operate, broadly 
speaking, as follows: Given a set of spectral energy distri- 
bution (SED) templates, the SEDs are projected onto pass- 
bands shifted on a grid of redshifts, from which the pre- 
dicted photometry, and hence colour, for each case (tem- 
plate type and redshift) is calculated. The observed colours 
are then compared against the model magnitudes, with a x^- 
minimisation over all bands used to determine the best-fit 
redshift and spectral type. The choice of templates is crucial 
for obtaining quality photo- 2;, so we will begin by describing 
our choice of templates. 



^^ For blue galaxies, [OIII] and Ha lines are both off the spec- 
trum (particularly in the vignetted corners where the spectral 
coverage is reduced), and for red galaxies, there is difficulty in 
confirming a redshift when there are no strong spectral features 
available, but only the forest of absorption lines. Both of these 
problems turn out to occur at similar redshifts in the DEEP2 
EGS survey. 
^^ http : //skyserver ■ sdss3 ■ org /dr8/en/ | 
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Figure 9. The six SED templates used in the ZEBRA photomet- 
ric redshifts. Four te mplates (labeled Ell, S bc, Scd and Im) are 
observed spectra from lColeman et al.l 198CI). and two (SB 2 and 
SB3) are synthetic starburst spectra 1 Kinney et al.lll996l ). Five 
templates were interpolated between the existing 6 templates, for 
a total of 31 templates which were used for ZEBRA photo-z es- 
timation. 



4.1 Templates 

Templates can be obtained from observation, or generated 
from stellar population synthesis (SPS) mo dels. We choos e 
the set of six SED templ ates as described in Bcnitoz (2000|), 
and used in ZEBRA (JFeldmann et al.l [2OO& ). They are 
shown in Figure [9l T he set consists of four templates from 
I Coleman etaP (igSO", CWW hereafter), which are SEDs ob- 
served across a wavelength baseline of 1400A < A < lOOOOA 
in the local (2 ~ 0) universe; they are supplemented by two 
synthetic blue spectra, us ing the GISSEL synthetic models 
(jBruzual fc CharlotI Il993h . While each of these templates 
have been linearly extrapolated into the UV and NIR wave- 
lengths based on different synthetic model^fj, these exten- 
sions are largely irrelevant to the optical filter set that we 
use (3OOOA < A < lOOOOA) for the redshift range relevant 
for our sample (z < 1.5). The CWW templates have been 
widely employed and are known to generate reliable pho- 
tometric redshift estimates, although trai ning and tweaking 
the templates improve s redshift estimates (jllbert et al.ll200(3l : 
iFeldmann et al.ll2006l). 



In particular, lllbert et al.l (|2006l ) find that these tem- 
plates work reasonably well for VVDS data, with a photo-2 
scatter of (Taz/{i+z) ~ 0.029 given deep photometry covering 
a similar range of wavelength to ours (u*g'r'i'z' and BVRI, 
with additional JK bands for 13 per cent of galaxies). With 
SDSS photometry, we expect a larger scatter because of the 
shallower photometry and lack of extra bands CJ 



http : //cas j obs . sdss . org/dr7/en/ 1 



http : //www . sdss . org/dr7/products/value_added/ Index. html | 



^"^ Th e SED extensions used here differ from those of lllbert et al.l 
1120061 ) , who use the same tem plate sets, but whose exten sions are 
updated with the more recent iBruzual &: Charloj 1120031 ). 
^^ We tested other template sets , such as the Poggiani tem- 
plates used bv lAbdalla et al.l 1I2OO 81) and PoUetta templates gen- 
erated by GRASIL jllbert et al.l2009.) . We find the CWW -I- Kin- 
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4.2 ZEBRA options 

For the difTerent template-fitting photo-z codes, their photo- 
z accuracy are very similar if the same set of templates is 
used, since most methods u se essentially similar bas ic pro- 
cedures of x^ minimization (jHildebrandt et al.ll201Cl ). How- 
ever, in addition to template choice, ZEBRA offers several 
options that enhances the photo-z calculations. Priors on 
the redshift distribution or other conditions may be ap- 
plied when calculating y^ . Additional features, such as it- 
erative photometry self-calibration, template optimisation, 
fe-correction tables, and different modes of photo-z and tem- 
plate type determination, are also available. We detail our 
choices of such ZEBRA options here. 



J^.2.1 ZEBRA parameters 

We ran ZEBRA in the maximum-likelihood (ML) mode. 
This was chosen over the Bayesian mode (BP), which allows 
for a Bayesian prior on the redshift distribution, per tem- 
plate type and redshift range. Although the latter method 
removes the photo-z bias to some extent, we find that it sig- 
nificantly adds to the scatter in photo-z, which is the quan- 
tity of most concern for our application. Because the original 
CWW templates arc known to generate photo-z's with bias 
(Fig. 5 of Fcldmann ct al. 2006), applying a "correct" red- 
shift distribution prior (the estimated redshift distribution 
of the given sample) will decrease the bias while increasing 
the scatter. 

We use the following set of parameters to obtain our 
photo-z's. The SED templates were interpolated with 5 new 
templates between each of the 6 given templates, resulting in 
a total of 31 templates. The redshifts are allowed to vary in 
steps of 0.0025 from to 1.5, with no smoothing of the filter 
bands. We have applied a prior of photo-z < 1.5, which is 
reasonable for single-epoch SDSS photometry (most galaxies 
within our magnitude limit are at z < 1.0). Si milarly, we do 
not c orrect for the Lyman-a IGM absorption (jMadaull 19951 : 
lMeiksin.2006. ) since it is only relevant for objects at z > 3. 



4-.2.2 ZEBRA self-calibration 

In principle, the template-fitting may benefit from use of 
a spec-z training set to calibrate the zero-point magni- 
tude of the photometry; here, this procedure is unnecessary 
as we find no offset for SDSS. This is as expected given 
the e xcellent photometric calibration (jPadmanabhan et al.l 
120081 ). The SED templates can be modified based on a spec- 
z calibration set, but we also find this template modifi- 
cation procedure unnecessary. While template calibration 



ney templates to be sufficient for the SDSS photometric samples. 
lAbdalla et al.l (2008) find that the Poggiani template set works 
well with the LRG samples. We confirmed their results for the 
redshift range 0.4 < 2 < 0.8 and the galaxy types in their sam- 
ple; however, it was not as successful as our fiducial choice out- 
side of that redshift range and colour selection. The Poletta tem- 
plates were required for good photo-2; estimates in lllbert et al.l 
1I2OO9I ) due to their use of IR photometry; in our case, with only 
optical data, we find no significant difference from the default 
CWW-I- Kinney templates. 



hav e been shown to work with samples of b etter photome- 
try (jllbert et al.ll2006l : iFeldmann et al.lbood '). our attempts 
at template modification suggested that we were only fit- 
ting to the photometric noise, as an improvement to the 
templates and photo-z error based on one calibration set 
did not translate to less scatter in the photo-z computed 
with other calibration sets. Moreover, the modifications to 
the templates did not look remotely physical. The differ- 
ence between our find ing and the improvements found in 
IFeldmann et al.l (|2006l ) most likely originates from the low 
signal-to-noise of our photometry. 



4.3 fc-correction 

ZEBRA, as a template fitting method, allows estimates of 
absolute magnitudes and stellar masses. To obtain these, 
we have utilised the fc-correction table (a list of theoretical 
magnitudes for each SED template at all redshifts; Eq. ([7]) 
below) generated by ZEBRA. This table is then modified 
according to the procedures below, to allow versatility in 
our fc-corrections. Factors of logjQ(l -|- z) that appear with 
different signs and/or prefactors in previous papers on the 
topic of fc-correction are elucidated here. 

fc-correction encapsulates the magnitude corrections be- 
yond the simple distance modulus, such that 



Mq =mR- DM(zgai) - /(-(zgai, T, R, Q). 



(5) 



Here, Mq is the rest-frame absolute magnitude in filter Q, 
rriR is the observed magnitude in band R, DM{z^a\) — 
51ogj^Q(DL/10pc) is the distance modulus (where Dl — 
DL{zgai) is the luminosity distance in parsecs), and 
K{zga_i,T,Q,R) is the fc-correction factor. The fc-correction 
depends on the galaxy redshift Zgai, template type T, the 
band R in which the magnitude is observed, and the desired 
band Q for the absolute magnitude. To minimise the degree 
of fc-correction in a sample, Q is often chosen such that it is 
the R band filter redshifted to the sample median redshift zq: 
i.e., Q ^^°R, or Q(A(1 -I- ^o)~^) = R{X). Then fc- corrcctions 
are calculated as (JHogg et al.ll2002l : iBlanton fc Roweis 20071 ) 



K{zg,u T, R, Q) 



2.51ogio(l + ^gal) (6) 

+ m(zgai, T, R) - m{z = 0, T, Q) 



where ZEBRA tabulates the AB magnitudes m for an array 
of redshifts Zgai : 



m(zgai, T, i?) = -2.51ogi 



/dAAFf[A(l-hZgai)-i]i?(A) 



JdX\g^^{X)R{X) 



(7) 

Here the quantity in square brackets is the flux ratio that 
defines the AB magnitude. Since the observed SED output 
of a galaxy of type T at redshift Zgai is 



(1 + Zgal)"'Fj[A(l -h Zgai)- V^ 



(8) 



where F^{X)dX is the rest-frame SED, we see that the 
2.51ogjo(l + Zgai) factor which appears in the first term 
of Eq. Q is from correction to the tabulated magnitude 
m(zgai) for the loss in photon energy (fiux) due to redshift. 
The reference SED relative to which AB magnitudes are 
defined is gt^{X), and [A gx'^iX)] ex A-\ since g^'^ is pro- 
portional to A"^. The arbitrary normalisations in both g ^ 
and F\ cancel out when we take the difference of m. 
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If Q is R' redshifted by zo, such that R' is one of the 
existing bands (i.e., Q{X{l + zo)~^) = R'{\)), then the tabu- 
lated values can also be used to obtain rh{z = 0,T,Q =^'H'): 



m(0, T, ^"i?') 



-2.51ogi 



-2.51ogio 



/dAAfJ(A)g(A) 
/dA A-i (5(A) 

/dA'A'i^f[A'(l + zo)-']i?'(A') 



(l + 2o)2/rfA' A'-ii?'(A') 

= m{zo, T, R!) + 51ogio(l + zo) 

where A' = A(l + zq). The fc-correction from the ZEBRA 
tabulated values are then 



J<(2gal, r, i?, ""i?': 



= 2.51ogi„(l + Zgai) - 51ogio(l + zo) 
+ m(2gai, T, -R)-m(zo, T, 7?(^0) 

where 7? and i?' can be any of the available filter bands. For 
the special case where R — R' and Zgai = zo, we have 



K{zo, T, R, ""i?) = -2.51ogio(l + zo) 



(11) 



which does not depend on the template type T by construc- 
tion, but is non-zero. Since the SED type T is uncertain in all 
photo-z's, the observed band R should be chosen such that it 
maximally overlaps with '°R! to minimize the T-dependence 
of the fc-correction. 



5 PHOTO-Z ACCURACY 

In this section, we report the ZEBRA photo-z accuracy on 
the SDSS DR8 photometry, using all of the spectroscopic 
calibration subsamples together. Since most of our galaxies 
are near the magnitude limit where the photometry is noisy, 
the derived photo-z's have a relatively large uncertainty. We 
consider dependence of photo-z scatter on template, and use 
this information to apply photo-z quality cuts. We have not 
found a significant correlation of the x^ values to the actual 
photo-z errors. 



5.1 Photo-z errors by template types 

Figure[TO]shows the photometric redshift error for the source 
catalog, for each template type separately. The errors are 
displayed as a function of photo-z, i.e., the observable; this 
format also clearly displays the different features in photo-z 
errors by template types. The dispersion around the median 
bias is quantified as 



2 



= {[Az - median(Az)]^) 



(12) 



and each panel shows ctaz for that template type. The photo- 
z failureqfj constitute ~2 per cent of the objects, which are 
not shown. The relative fraction of template types in the 
photo-z sample show that approximately a quarter each are 
of the Ell, Sbc, and Scd galaxy types; about 10 per cent 
are of the Im-type, and <10 per cent are classified as either 
SB2 or SB3. The lens sample shows a similar distribution of 
galaxy types. 

The galaxies classified as starburst (SB2 and SB3) SED 




^^ Photo-2's at the ZEBRA prior limits, Zphot 
considered photo-z failures. 



and 1.5, are 



Figure 10. Photometric redshift errors in the source catalogue 
by template types (EII, Sbc, Scd, Im, SB3, and SB2), as a func- 
tion of the photometric redshift (the observable) . The number in 
parenthesis indicates object counts. The median (bold solid line) 
and the 16 and 84 percentile range (bold dashed line) of Az is also 
shown. The photo-z scatter is given as the standard deviation of 
Az from the median. Only 50 per cent of the points are shown 
for clarity. We impose a photo-z quality cut by removing all SB2 
and SB3 objects from our bias analysis, with < 10 per cent loss 
of objects. 



types show the largest photo-z errors. Our photo-z quality 
cut discards the SB2 and SB3 template type objects along 
with the ~2 per cent photo-z failures. We lose <10 per cent 
of our sample with this cut. 

For the remaining four template types, the width of the 
photo-z error (dashed lines, Fig. llOp in general becomes nar- 
rower as the galaxy becomes brighter, down to ctaz ~ 0.09 
for r < 20.5. In general, it is believed that red galaxies (Ell) 
have better photometric redshift accuracies compared to the 
bluer galaxy types, due to the high contrast below and above 
the (rest-frame) 4000A break. However, this statement is 
only true when there is a strong photometric detection in 
the bands above and below the break. The large photo-z 
errors in the Ell and Sbc samples at Zphot < 0.2 are due to 
noisy photometry; this broad dispersion stems from the fact 
that photo-z depends on the colour difference between the 
g and the ri bands at low Zphot to probe the 4000A break; 
but at dim ri magnitudes, the g magnitudes for the red Ell 
and Sbc galaxies are even dimmer, resulting in a photo-z fit 
that is highly sensitive to noise. 

Figure [11] shows results for the full samples (all tem- 
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Figure 11. Photo-z as a function of spectroscopic redsiiifts for 
the SDSS lens (top) and source (bottom) catalogues, after photo- 
z quality cuts have been applied. The median (bold solid line) and 
the 16 and 84 percentile range (bold dashed line) of A.z is also 
shown. The error uncertainty is given as the standard deviation 
of l\z/{l + Zspec). Only 33 per cent of the points are showrn for 
clarity. 
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Figure 12. SED template type uncertainty in the source cata- 
logue (the lens catalogue show similar scatter). The ordinate is 
the "true" template ID, while the abscissa is the best template 
type estimate derived from the photo-2 run, and is our observ- 
able. The thick solid line and thick dashed lines show the median 
and the upper/lower 68 percentile scatter, respectively. The 6 ba- 
sic CWW-|-Kinney template types correspond to template type 
IDs 0, 6, 12, . . . 30, with 5 interpolated SED templates in be- 
tween. Our SED type classification bins are delineated by the 
straight, dotted lines. (The points are shown with random frac- 
tions of S [—0.5,0.5) added to better illustrate the density; the 
template type IDs are all integers.) 



particular. Figure [TT] shows bimodality in the photo-z distri- 
bution for 0.2 < Zspac ;$ 0.4. This feature is seen in all tem- 
plate types, and is an artifact of the gap between the g and 
r-bands at 5440 A, or equivalently, when the 4000 A break 
transitions from g into r at z ~ 0.36. The gap between the g 
and r bands causes uncertainty in the photometric redshift 
when Zapec = 0.36 even with perfect photometry; with high 
photometric noise, the error in the photo-z becomes bimodal 
centered at this redshift. 



plates) using the more standard format for comparison with 
other photo-z literature: the photo-z's are shown as a func- 
tion of the true redshift, and scatter is quantified as 



0"Az/(i-i-z) = var[Az/(l + Zspc 



(13) 



where Az ; 



Zphot 



The photo-z error is large, where 



''"Az/(i-i-z) = 0.096 and 0.113 for the lens and the source sam- 
ple, respectively, averaged over the whole sa mple. In compar- 
ison, ZEBRA template-fitting photo-z's in iFeldmann et al.l 
l|2006t ) achieve accuracies of o-az/(i+z) < 0.03 when applied 
to Iab < 22.5 COSMOS galaxies with u* BVg'r'i'z'Ks pho- 
tometry, thus suggesting that the systematic fioor due to 
limitations of the code and template set are low, and our 
errorbars are dominated by photometric noise. We find that 
the scatter is a function of magnitude, and degrades rapidly 
for r > 21. For our source catalogue, however, the degra- 
dation at the faint end is not so severe (crAz/(i+z) ~ 0.13 
for 21 < r < 21.8); this implies that the better resolution 
required for the source sample yields more reliable photom- 
etry. 

The photo-z scatter is asymmetric about the median. In 



5.2 Uncertainty in template types 

It is important to understand the error in template type des- 
ignation, not only because we define our photo-z quality cut 
based upon template type, but also because we want to un- 
derstand and minimise the errors in fe-correction (Sec. I4.3|l 
and stellar mass estimates (Sec. l6.2TT|) . Figure [121 shows the 
template designation errors for the source catalogue. The 
"true" template types were determined by asking ZEBRA 
to fit to the best template t ype , given the known Zspec, and 
is plotted on the vertical axi^fj; the horizontal axis is the ob- 
servable (the template type derived from photometry alone, 
when determining the photo-z). There are 5 interpolated 
SED templates in between the 6 basic CWW-|-Kinney tem- 



^^ The "true" template types in this section are true to the extent 
that the ZEBRA template set provides an accurate way to clas- 
sify galaxy templates. Thus it is meant to represent truth in the 
absence of what photometric noise does to the ability to estimate 
a redshift. 
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plate types; our SED type classification bins are delineated 
by the straight, dotted lines. 

The median of the scatter lies on the slope = 1 line, 
and the scatter width from the median in the range of 
Ell/Sbc/Scd/Im SED types is ± 3 units, which is approxi- 
mately the half-width of the template type bins. Hence we 
see that the template type error is reasonably small for the 
chosen template categorisation. In contrast, very few ob- 
jects fit to the interpolation between Im to SB2 or SB2 to 
SB3 types, with large uncertainty in the template types in 
these bins. The lens catalogue template type uncertainty 
show similar scatter. 

However, even when using the spectroscopic redshifts, 
it is possible to have the wrong template, either because the 
template set may not describe the true galaxy SED, or due 
to the relatively large noise in the magnitudes which may 
cause the selection of the wrong template even when given 
the true redshift. This additional template type uncertainty 
is not accounted for in Fig. 1121 



6 BIAS IN LUMINOSITY AND STELLAR 

MASS 

To obtain a weak lensing signal with sufficient signal to 
noise, lens galaxies with similar properties (such as lumi- 
nosity and stellar mass) are stacked. Absolute magnitudes 
and stellar mass, used to assign lens galaxies to each stack, 
inherit scatter from photo-z errors. With our lens calibra- 
tion set we can measure this scatter, and use the width of 
the scatter to determine the binning resolution for stack- 
ing. Moreover we can test for biases in the absolute mag- 
nitudes and stellar masses, and figure out how to shift our 
bins around to account for this type of error. 



6.1 Absolute magnitude 

The bias in absolute luminosity due to photo-z errors is a 
combination of two effects: the error in the luminosity dis- 
tance, and in the fc-correction (Sec. 14. Sp . All magnitudes 
are fc-corrected according to SED template types. We define 
our absolute magnitude band as the r magnitude redshifted 



to z 



0.3 C'^Mr), so that in Eq. ([T0|, 7? = R' 



and 



zo — 0.3, and where the observed r magnitude was chosen 
for its high S/N. These choices minimise the degree of k- 
correction for our lens sample, and most of the correction 
is dominated by the distance modulus, as demonstrated in 
Figure 031 Here, the theoretical absolute magnitude errors in 
"'^Mr (distance modulus -I- fc-correction) are shown in solid 
curves, for objects at Zspec= 0.1, 0.2, . . . , 0.6, as a function 
of Zphot. The dashed lines show the magnitude modifica- 
tion from the distance modulus component alone; the fact 
that they nearly coincide with the solid lines indicates that 
the distance modulus error is more important than the un- 
certainty in fc-corrections. The points show the distribution 
of galaxies from our photometric lens catalogue (r < 21). 
Although the fc-correction is in general small, it becomes 
somewhat more significant for Ell objects at high Zphot or 

Zspec . 

The size of the scatter in "'^Mr dictates the appropriate 
luminosity binning width. We further bin the lens objects 
by redshift bins of width Azphot ~ 0.1, up to Zphot ~ 0.6, to 




0.2 0.3 
photo-i 



Figure 13. The theoretical "'^A'/r absolute magnitude error 
curves (distance modulus -I- fc-correction) for objects at Zspec= 
0.1,0.2, . . . , 0.6 (solid curves, from bottom to top), as a function 
of Zphoti for each template types (Ell, Sbc, Scd, Im). The dashed 
line is the distance modulus component alone. If the photo-z is 
below the spec-z (top half of each panel), the estimated abso- 
lute magnitude is too large (fainter) than the real magnitude. 
The correction is more severe at smaller photo-z's. The points 
show distribution of objects from our photometric lens catalogue 
(r < 21). 
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Figure 14. Absolute magnitude bias due to redshift errors, in 
Zphot '^'ud ^'^Mr bins, for the Ell SED type. The 6 panels cor- 
respond to different Zp^ot bins; the points show individual ob- 
jects, while the thick solid and dashed lines are the median and 
16/84 percentile of the points in a given ^'^Mr bin. The ver- 
tical dotted lines correspond to the true '^■^Mr of the median 
Ap-'^A/r], and these lines are labeled with the true "'^Afr values. 
The large bias and scatter at low photo-z (zpijot< 0.2), low lumi- 
nosity C'^Afr > —21) bins are due to the large photo-z errors in 
those bins for this SED type. The numbers in parenthesis indicate 
the galaxy counts in the panel. 
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Figure 15. Absolute magnitude bias due to redshift errors, in 
Zpliot and "■'^Mr- bins, for tiie Sbc SED type (see Fig. 1141 caption 
for description). 



Figure 17. Absolute magnitude bias due to redshift errors, in 
•2phot a-nd ^■'•^Mr bins, for the Im SED type (sec Fig. 1141 caption 
for description). 
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Figure 16. Absolute magnitude bias due to redshift errors, in 
Zphot and ^'^Mr bins, for the Scd SED type (see Fig. 1141 caption 
for description). 



determine the binning width as a function of ^phot and SED 
type. The actual distribution of absolute magnitude error 



A["^^M,] 



''Mr(zphot) 



Mr{ztT\xa) is determined by 



the photo-2 distribution and error distribution of the photo- 
metric lens sample. Figures [T4l through 117! show the median 
and scatter in the error A['^''^i\fr] for the four SED template 
types T that we use for our analysis, as a function of the es- 
timated absolute magnitude ''■^Mr(2phot) (the observable), 
where sub-panels in each figure are binned by the photo- 
metric redshifts. 

The smallest scatter is seen in the Ell type galaxies 
at high redshift (z > 0.3) and bright absolute magnitude 
{^'^Mr(z-p\^o\) < —20). This is the binning range most rele- 
vant to current SDSS LRG studies. These bright, red galax- 
ies show the smallest A[''^^Mr] scatter of om^ ~ 0.7 mag- 
nitude from the median. The large scatter in the low Zphot, 



low luminosity region, A[°^^Mr.] > 2, comes from the broad 
redshift uncertainty in the red SED types at low r magni- 
tude, as discussed in Sec. 15.11 Unlike LRGs, dim red objects 
with "■^Mr(zphot) ~°'^Mr(2spcc) > -21 do not have reliable 
magnitudes, because the photo-2 errors are worse (Sec. 15.1] ). 
For other SED types and photo-z bins, a typical magnitude 
scatter is about om^ ~1: and so the bin size cannot be much 
smaller than this. We also note that where the calibration 
sample distribution is sparse, the binning cannot be reliable. 
For this reason, the redshift bin < Zphot < 0.1 is probably 
best omitted, for lenses of all SED types. 

The vertical dotted lines in Figures I14H17I show where 
stacking the lens galaxies by absolute magnitude binning 
would be meaningful where the average cro.ajv/^ (averaged 
over the photo-z bin) are shown in each panel. The brightest 
bins ("'^Mr < —21) are accessible to the Ell and Sbc types at 
high Zphot, where the faint end ("'^Mr > —19) are accessible 
to Scd and Im objects in the 0.1 < Zphot < 0.2 bin. 

An additional consideration in comparing magnitudes 
across different redshifts is the luminosity evolution. Be- 
tween the lowest and highest redshift bin, evolution ac- 
counts for a difference in ~0.5 magnitudes. This is esti- 
mated according to the luminosity function (LF) evolu- 
tion, based on several LF studies in a similar redshift range 



dWolf et al.ll2003l:lGiallongo et alllioosl : I Willmer et al.ll2006l : 
iBrown et al.l l2007l : iFaber et all \2Qm . We adopt Q 



^ . 1.2 

mag/redshift for all templates, which appears to be consis- 
tent with all of these studies, and where Q is the redshift 
evolution slope as defined in 



M,(zo) = M,{z) + {z- zo)Q 



(14) 



where Af« is the magnitude corresponding to L* in the lumi- 
nosity function, z is the object redshift, and zq is the stan- 
dard redshift to which we normalise our magnitude binning. 
These studies are based on B-band luminosities, but °\ and 
B are close enough (they partially overlap at 4300-4700A); 
hence the correction from the colour evolution is expected 
to be small. If we make the simplifying assumption that 
the luminosity rank ordering does not mix as the luminosity 
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function evolves, then we can define magnitude bins that 
shift according to the slope Q as a rough approximation to 
the luminosity evolution. This correction is a small effect 
compared to the shifts induced by photo-z error bias. 

6.2 Stellar mass 

6.2.1 Estimating stellar mass 

Stellar mass can be estimated from the combination of (1) 
rest-frame colours and (2) absolute mag nitude (di s cussed 
above in Sec. 16. ip . using correlations by iBell et al.l (|2003l ) 
(via an intermediate step of estimating a stellar mass-to- 
light ratio, M,/L). The colou r — Mt,/L r elation is defined 
using the rest frame colour in iBell et al.l (2003); and SED 
types are used to convert colours and magnitudes into the 
rest frame. As excessive fc-corrections are susceptible to er- 
rors in photo-2 and SED types, we choose the optimal ob- 
served band which minimises the fc-correction to the rest 
frame band. We fc-correct the observed i magnitude into the 
"■"V band rest-frame absolute magnitude and the observed 
r — i colour to °°{g — r) (Sec. l4.3|l . The °°r band rest-frame 
is chosen since its wavelength range approximately corre- 
sponds to that of the observed i band at the lens median 
redshift of zo = 0.3, which is the reddest band available in 
SDSS with high photometric S/N (red luminosities require 
less correction to the M,/L ratio). Similarly the °°((? — r) 
colour is chosen for the higher S/N in the observed r and 
i bands. Then, t he M,/L ratio is obtained from Table 7 of 
iBell et al.l (|2003l ). 



logio(M./L,) = -0.306 + 1.097 ^ \g - t)\ 



(15) 



where M*/Lr (stellar mass M, to "'"r-band luminosity ra- 
tio) is in solar units. The r band luminosity in solar units 
is 



logio ^r = -0.4 ("■ Vr - 4.64) 



(16) 



where {Mr )(>, — 4.64 is the r-band absolute AB magnitude 
of the Sun IjBlanton fc RoweiauOOTI ). Hence the stellar mass 
in solar units is 



logio A'/. = 1.550 + 1.097 f'°{g - r)] 



0.4 \°-°Mr] 



(17) 



These numbers quoted above are based on a diet Salpeter 
IMF; for future science work we can modify for this by rescal- 
ing all stellar masses by a constant factor. 

We now consider biases in these stellar masses. We start 
with biases in the intermediate step of constructing the stel- 
lar mass-to-light ratio. 
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Figure 18. The scatter in the stellar mass-to-light ratio (M«/L) 
with respect to rest-frame r-band luminosity in solar units, as a 
function of photo-2, for the different SED types, determined from 
photo-2 fitting. The thick horizontal line shows the theoretical 
Mt/ L based on the rest- frame colour, and is constant for a given 
SED type. The points are the M«/L estimated from the "true" 
template type (best-fit template given the spec-z) fe-corrected us- 
ing the true redshift (zgpec), and the solid/dashed lines show the 
median and 68 percentile scatter, respectively. The Af, /L deter- 
mined from Zpiiot has intrinsic scatter of ~0.1 dex (not shown). 
The scatter in Mt/L is small compared to the bias in the lumi- 
nosity L. 



"true" template typqfj and fc-correct using the true red- 
shift. The fc-corrected colour ""((?- r), and hence the derived 
Mf/L, then show some bias and scatter, as seen in Figure fTSl 
The constant (as a function of Zphot) predicted M,/ L for the 
given template type is shown for comparison. The points 
show the distribution of the lens catalogue objects, and the 
solid/dashed lines show the median and 68 percentile scat- 
ter, respectively. There is still little variation in the M^/L 
evident, because the rest-frame colours are identical for a 
given template type, and the "true" template type deviates 
little from the estimated one. The magnitude of the scat- 
ter shown is similar to the uncertaint y in the M,/ L re lation 
(Eq. [151) itself, which is ~ 0.1 dex (|Bell et al.| [ 2003l ). The 
true errors in Mt,/L might be expected to be larger than 
that shown here, given that the actual galaxy SED is un- 
known. However, we expect the error in the luminosity to 
dominate, discussed below. 



6.2.2 Bias in stellar mass-to-light ratio 

In this section, we consider errors in the inferred M-t/ L us- 
ing photo-z, starting from the basic assumption that these 
are signifi c antly larger than intrinsic uncertainties in the 
iBell et al.l l|2003r ) prescription for estimating stellar mass 
when using spectroscopic redshifts. This statement ignores 
underlying uncertainties in the stellar IMF, which plague all 
stellar mass estimates. 

There is very little bias and scatter when we plot the 
Mt/L obtained from the rest-frame ""(ff- J') colour by SED 
template types. This is because for a given SED template, 
the rest-frame colour is a constant. Instead, we take the 



6.2.3 Bias m stellar mass 

While the scatter in Mt/L is of order 0.1 ~ 0.2 dex, the 
scatter in the luminosity is ~ 0.4 dex (converted from ~1 
mag scatter, see Eq. ll6p . mostly dominated by the distance 
modulus (Sec. 16.1"]) . Figures [1914221 show the bias and scat- 
ter in stellar mass estimate in the four template types, as 
a function of estimated stellar mass, for the 6 photo-z bins. 



^■^ The "true" template type is the best-fit SED template given 
2spec, which is unbiased but scattered with respect to the observed 



template type; see Fig. 1121 
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Figure 19. Stellar mass errors Alogjo[A^*/A^0] = 

Iogio[M*/Af0](zphot) - logi(,[M./M0](2,pec) in the Ell SED 
type, as a function of the estimated stellar mass based on Zphot 
(the observable) . The distribution of the lens catalogue are plot- 
ted as points, where the median and 68 percentile scatter are 
shown as the thick solid and dashed lines, respectively. Each panel 
correspond to different Zphot bins. The top axis of the top panels 
show the nominal stellar mass bin of 0.3 dex, or where the bin 
is twice as massive as the previous one. The vertical dotted lines 
show the bias-corrected binning. 
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Figure 20. Stellar mass errors in the Sbc SED type. 



The errors (vertical cixis) are plotted as a function of esti- 
mated stellar magnitudes \og-^Q[Mt/MQ]{zphot) (the observ- 
able), where the nominal stellar mass bins are in widths of 
0.3 dex (i.e, each bin is twice as massive as the previous bin). 
These bins are labeled by the index 1 through 8, as seen on 
the topmost panels in each Figure. The vertical dotted line, 
labeled by stellar mass bin numbers, are the corrected stellar 
masses based on the median values. 

We see the general trend in the absolute luminosity 
bias and scatter repeated in these plots, e.g., the Ell types 
show small scatter at high stellar mass and high Zphot, and 
the lowest Zphot bin (0 < Zphot < 0.1) is too sparse. The 



Figure 21. Stellar mass errors in the Scd SED type. 
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Figure 22. Stellar mass errors in the Im SED type. 

logjfl[A/*/MQ] scatter criogM, , which is defined with respect 
to the median, is shown in each plot, and indicates small 
scatter for most bins. In fact, most bins show scatter smaller 

than the nominal 0.3 dex stellar mass bin. 

W e note that previous studies (e.g., iDutton et al.l 
l|201Cl )) have compared the photometry-based stellar mass 
estimate employed here with spectra-based estimates. On 
average, the two agree for higher stellar masses, but deviate 
systematically at lower stellar mass. Such systematic devia- 
tion can be corrected accordingly with a stellar-mass depen- 
dent correction. 



7 BIAS IN THE GRAVITATIONAL LENSING 
SIGNAL 

With these photometric redshifts, their scatter and their ef- 
fects on stellar mass and absolute magnitudes quantified, 
we can turn to the bias and scatter induced in the gravi- 
tational lensing signal due to photo-z use for sources and 
lenses. As the conversion of lensing shear to mass depends 
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upon the redshifts of the lens and source (see Eq. [T|, errors 
in these redshifts propagate into errors in the weak lensing 
measurements. We consider three cases, (a) the lens redshift 
is known and photo-z is only used for the source, (b) the 
source redshift is known and photo-z is only used for the 
lens, and (c) photo-z's are used for both lens and source. 
We deal with these in turn below, focusing on biases in the 
lensing signal calibration; in principle, there is also a blur- 
ring of information in the transverse direction once photo-z 
are used for lenses; however, a typical 10 per cent photo- 
z error does not have a significant effect on the shape of 
the lensing profiles, so we neglect this effect here. We be- 
gin by reviewing and then expanding upon the methods by 
iMandelbaum et al.l (|2008l ) to calculate lensing bias. 



7.1 Methods 

7.1.1 Bias 

In the absence of photo-z errors, the observed lensing tan- 
gential shear 74 can be related to the lens surface density 
contrast AE via 



AE = Se7t 



(18) 



where the proportionality constant, the critical surface den- 
sity, was defined in Eq. ((T}. 

The redshift calibration bias is the misestimation of 
AE due to the photo-z scatter. To estimate this bias, we 
first consider the method of estimating AE, via a weight- 
average over the tangential component of the measured 
source galaxy shapes 7t , 



AE: 



Ei ^1 



(19) 



Here, the summation is over j lens-source pairs, E^j is the 
critical mass density for the jth lens-source pair, and the 
tilde indicates estimated values (using photo-z's where ap- 
plicable). The optimal weight is 



(jMandelbaum et al.l I2OO8I ) . The analytic curve is fit sepa- 
rately for every calibration subsample that is used. 

The bias (Eq. l21|l is for a single lens redshift. If we want 
to estimate the average bias from using photo-z for Zsrc over 
all lenses with known zicns, we need to average over the lens 
redshift distribution, including source weight factors. 



{b.) = 



/dZlensP(2lens)uii(ziens)6z(2lens) 
/ dZlons p(zicns) W! (ziens) 



The weight for a given lens redshift w;(zions) is 

W'i(2lons) = -D^^(l + Zlcns)"^ 5Z "'*' ' 

k 



(22) 



(23) 



Here the summation is over all source galaxies k estimated 
to be beyond the lens redshift zicns, and zions may be spec- 
troscopic or photometric. The reason for the prefactors be- 
fore the summation is that our method of estimating these 
sums using a calibration sample of fixed area does not really 
correspond to how lensing signals are actually measured. 
Typically they are estimated within some fixed physical or 
comoving aperture, which means that lenses at lower red- 
shift will use sources from an effectively larger area on the 
sky, and therefore get greater weight according to the square 
of that angular diameter distance to the lens redshift. The 



factor of (1 4- zio 



assumes the use of a fixed comoving 



aperture. Note that this effe ct was incorrectly neglec ted in 
the previous analysis (eq. 6 of IMandelbaum et al.l2008r ). For- 
tunately, it is only significant when averaging over a broad 
lens redshift distributio n . In s ubsequent papers relying on 
the IMandelbaum et al. 1120081) lensi ng signal calib r ations , 
including [Mandelbaum et al.l (I2OO8I'): iReves et al.l (|2008| ): 
IMandelbaum et al.l (|2009l . 120101 ): [Schulz et al.l (|2010l '). the 
change in the lensing signal calibration when including this 
additional redshift weighting factor is typically 2 per cent, 
which is approximately the size of the quoted systematic 
uncertainty in the redshift calibration, and typically 20 per 
cent of the statistical error. 



E2 . 



{e 



2 -f 0-2 

rms ' ^ e ' 



(20) 



The quantities added in quadrature are erms, the shape noise 
of the source ensemble, and Uej, the shape measurement er- 
ror of the jth galaxy (per single component of the shear). 
Assuming that the only calibration bias occurs through the 
use of photo-z via thejjuantity Ec, the redshift calibration 
bias 6z is the ratio of AE to the true signal AE. By substi- 



tuting 7^^ 



U) 



^SE-j, 



we find 



6z(zions) + l = -T^ 



AE 
AE 



E,^jS-JE, 

E,- wj 



(21) 



the weighted sum of the ratio of the estimated to true critical 
surface density. 

In the actual bias estimation, each galaxy j is further 
weighted to "smooth" out the LSS in the calibration sam- 
ple. To do so, we fit the redshift histogram to the analytic 
curve Eq. ((Sj, as seen in Figs. [8]and[7| the LSS weight for 
all galaxies in a particular histogram bin i is then the ratio 
of the number of galaxies according to the fit distribution. 



to the real number of galaxies, i.e., wlss 



N. 



(model) 



/N, 



7.1.2 Source-lens pair (statistical) incompleteness 

The bias as described above does not indicate the statisti- 
cal loss in the number of valid source-lens pairs, or in the 
deviation from optimal weighting due to the use of photo-z. 
Such contributions can be described in terms of the purity 
and completeness of the lens-source pairs, and the change 
in weights from using photo-z instead of the true redshift. 
Purity is the fraction of lens-source pairs that are truly 
lensed with all "valid" pairs based on photo-z, where each 
pair is weighted according to Eq. (|20|l . Completeness is the 
ratio of the sum of all weights of the lens-source pair en- 
semble, calculated based on photo-z to that based on the 
true redshift. The loss of statistics due to low purity and 
completeness translates to higher uncertainty (variance) in 
AE. The efficiency is defined as the ratio of the variances, 
where the smaller, ideal variance is attained when all red- 
shifts are known, and the estimated variance obtained from 
using photo-z's is always larger than the ideal variance due 
to the two previous issues (low purity and/or completeness) 
plus deviation of the weighting from the optimal weighting. 
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Figure 23. Lensing signal calibration bias from using photo-z for 
ZsTc, given a known lens redshift zions- The thin solid lines show 
bias in the EGS, VVDS, COSMOS, PRIMUS-02hr and PRIMUS- 
23hr fields, and the weighted average of the 5 solid curves is shown 
as the thick dashed line. Of the 5 fields, the COSMOS field bias 
has been empirically corrected (see text); the pre-correction bias 
is shown as the thin dot-dashed line. 



(24) 



This efSciency is given by (|Mandelbauni et al.ll2008l ') 

Ideal Var(AS) _ (^ ' 
Real Var(AE) ^ (^ w) (X^ '< 

where the weights are defined in Eq. (|20[1 . and the tilde 
indicates values estimated based on photo-z's. These will 
be illustrated for the case of both source and lens photo- 2: 
(3L2) below. 



7.1.3 Bias uncertainties from LSS 

The uncertainty in the calibration bias estimated using 
Eq. (|2ip is determined not only by the size of the calibration 
sample (Poisson statistics) but also by i ts LSS fluctuations. 
To est imate the size of this uncertainty, iMandelbaum et al.l 
1(2003) used a modified bootstrap resampling method to esti- 
mate the uncertainty in the true source redshift distribution, 
and hence in the bias curve, 6z(zicns), where Zieus is known. 
In this paper, we take a simpler approach. Since we have 
several calibration subsamples representing independent re- 
alizations of LSS, we obtain a bias curve 6z(zicns) for each 
subsample, take their average as the bias, and the standard 
deviation of the mean as the systematic uncertainty due to 
LSS (see, e.g. Fig. [23j. 

Given the formalism described in this section, we now 
turn to the case of source photo-z, lens photo-z, and both 
combined. 



7.2 Lensing bias in source photo-z 

First, we consider the case where the lens redshifts are 
known to high accuracy, but the source redshifts rely on 
much more uncertain photo-z. Figure [251 shows the lensing 
bias 6z(ziens) for using photo-z in the new SDSS source cat- 
alogue as a function of the known lens redshifts ziens. Bias 



for each the 5 source calibration subsamples (EGS, VVDS, 
COSMOS, PRIMUS-02hr and PRIMUS-23hr) are shown as 
thin solid lines. Of the 5 fields, the calibration bias in the 
COSMOS field has been empirically corrected to account 
for excessive sky noise, as will be described shortly. Prior to 
calculating the bias, all subsamples have the photo-z qual- 
ity cuts applied (Sec. I5.l"|) . and the "excess" low-resolution 
galaxies have been omitted in the VVDS, PRIMUS-02hr and 
23hr fields (Sec. 13.23]) . The bias 6z(zicna) can be as low as 
—0.4 at ziens ~ 0.5, but the uncertainty in the bias is less 
than 0.03; i.e., —0.40 ± 0.03. In general, source photo-z's 
with large scatter tend to pull the bias more negative. 

As noted previously, it was impossible to augment the 
COSMOS sample to correct for the loss of observable galax- 
ies due to excessively high sky noise (Sec. I3.2.ip : instead, 
it is necessary to derive an empirical correction to the bias. 
The redshift calibration bias curve for the COSMOS cal- 
ibration sample, before correction, can be seen as the dot- 
dashed line in Fig. 1231 clearly it is quite discrepant compared 
to the other samples. We have empirically compensated for 
this effect using the following procedure: We used the fact 
that our PRIMUS-02hr and 23hr calibration samples are 
on stripe 82, where there are many SDSS observing runs 
with different conditions. Thus, we chose runs overlapping 
those two regions with seeing that is comparable to in the 
COSMOS region, but with two different values of sky noise: 
typical, and 20 per cent larger (as in the COSMOS region). 
We carry out the bias calculations independently for each 
of the observing runs with different conditions, including 
matching against the spectroscopic training sample, impos- 
ing the lensing quality cuts, and obtaining photo-z's from 
the photometry. We calculate the resulting lensing signal 
bias 6z(zicns) from the photo-z's for the two runs separately, 
and then use the difference between to two to correct for 
the COSMOS sample bias. It is reassuring that after apply- 
ing this correction, the calibration bias for the COSMOS 
sample is reasonably consistent with that for the other four 
calibration samples. 

7.3 Lensing bias in lens photo-z 

Next we calculate the bias in the gravitational lensing sig- 
nal when the source redshift is known, but photo-z is used 
for the lens redshift. Unlike the case above (Sec. 17.2]) where 
zions is known but photo-z is used for source redshifts, such 
a scenario is not relevant for real observations. However, this 
exercise will allow us to study the general response of the 
lensing bias 6z to errors in lens photo-z, and may suggest 
certain scenarios in which using photo-z for ziena would be 
inadvisable. As this is just an exercise, we do not consider 
the purity, completeness and efficiency; only the lensing cal- 
ibration bias is shown. 

Figure [24] shows the bias 6z(zsrc) as a function of the 
known Zsrc- We initially found that all calibration subsam- 
ples agree well except for EGS; this deviation was due to the 
lack of objects in the z ~ 0.3 bin, as discussed in Sec. 13.41 
This discrepancy is easily remedied when the deficient bin 
in the EGS sample is replaced with, e.g., objects from the 
VVDS survey of the same redshift bin (both the original and 
corrected EGS biases are shown in the plot). Note that this 
correction was not required in Sec. 17.21 when using source 
photo-z. The observing conditions that affect the r < 21.8 
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Figure 24. Bias in the gravitational lensing signal due to photo- 
2 usage for lens redshifts, as a function of known source redshift, 
6z(2src). The different panels show Zphot bins of lens redshifts; 
the number in parenthesis indicates the number counts of galax- 
ies in the spectroscopic lens calibration sample in this photo- z 
bin. The thin lines are the biases in the six individual lens cal- 
ibration subsamples, and the thick line is the overall bias from 
the combination of all calibration subsamples. The thin dashed 
line indicates the EGS sample before the deficient z ^ 0.3 bin 
has been corrected for (Sec. I7.3|l . We discard the lowest redshift 
bin, < 2phot < O-li since the different calibration subsamples 
suggest inconsistent values for the calibration bias. 



source calibration sample are not relevant to the lens r < 21 
sample because (a) it does not have resolution cuts, and 
hence is not affected by the variation in the seeing, and (b) 
it has a brighter flux limit, such that the sky noise has min- 
imal effect on the number counts (Sec. l3.2T2|) . 

For all zions bins, we see that the lensing bias converges 
at high Zsrc to a reasonable value for most Zphot bins ex- 
cept the lowest lens redshift bin < 2phot < 0.1. We also 
see that the bias &z(zsrc) varies rapidly close to the lens 
redshift, and also shows the largest difference between the 
calibration subsamples (perhaps reflecting a sensitivity to 
LSS fluctuations in the lens number density or overall char- 
acteristics such as domination by red vs. blue galaxies). The 
characteristic width of this high-variation region is approxi- 
mately Azphot ~ 0.1, implying that a minimum lens-source 
separation of 0.1 in redshift can help stabilise the bias and 
minimize its associated systematic uncertainty. In reality, 
with only photo-z for both zions and Zsrc, such a clean cut 
is impossible, and there may not be significant benefit when 
applying a minimum separation cut. We explore this issue 
further in the next section. 
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Figure 25. Top left: Bias hz in the galaxy-galaxy lensing sig- 
nal due to photo-z usage for both the lens and source. The lens 
galaxies have been binned according to photo-2 bin width of 
A^phot = 0.1. All 5 calibration subsamples (EGS, VVDS, COS- 
MOS, PRIMUS-02hr and PRIMUS-23hr) have been plotted to 
show the LSS-induced scatter. Dots: no pair rejection applied. 
Crosses: Minimum photo-2 separation of 0.1 between lens and 
source required. Top right: Efficiency (degradation in variance), as 
defined in Sec. 17.1.21 Bottom left and right panels: Completeness 
and purity, respectively, as defined in Sec. 17.1.21 The weighting 
scheme keeps the lensing signal relatively pure, while the degra- 
dation in the variance (statistical loss) comes mostly from the 
incompleteness of the valid lens-source pairs due to the use of 
photo-2 errors. 



Table 3. Lensing signal bias and its uncertainty, for photo-2; lens 
and source redshifts, from the bias panel in Figure [25l 



Lens photo-2 bin 


Lensing bias 


(COSMOS noise bias) 


Using all samples 






0.05 


0.513 ± 0.055 


(-0.027) 


0.15 


-0.061 ± 0.025 


(-0.035) 


0.25 


-0.211 ± 0.027 


(-0.048) 


0.35 


-0.227 ± 0.029 


(-0.068) 


0.45 


-0.178 ± 0.034 


(-0.105) 


0.55 


-0.183 ± 0.069 


(-0.167) 


With minimum photo-2 separation 


of 0.1 


0.05 


0.521 ± 0.058 


(-0.027) 


0.15 


-0.063 ± 0.019 


(-0.041) 


0.25 


-0.219 ± 0.025 


(-0.049) 


0.35 


-0.244 ± 0.031 


(-0.074) 


0.45 


-0.220 ± 0.033 


(-0.115) 


0.55 


-0.249 ± 0.062 


(-0.167) 



7.4 Lensing bias, lens and source photo-z 

One of the main goals of this paper is to lay the ground- 
work for galaxy-galaxy lensing using photo-2: for sources and 
lenses. As mentioned in the introduction, this will dramati- 
cally increase the range of SDSS galaxies accessible for study 
with galaxy-galaxy lensing, both in the direction of lower lu- 
minosity at fixed redshift, and also going to higher redshift 
samples at fixed luminosity. Figure E51 shows the lensing sig- 
nal bias, efficiency, completeness and purity when photo- 2; 
lens-source pairs are used. The bias values have been tab- 
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ulated in Table [S] The lens galaxies have been binned by 
photo-2 only. All 5 calibration subsamples (EGS, VVDS, 
COSMOS, PRIMUS-02hr and 23hr) have been plotted to 
show the LSS-induced scatter. Each subsample has been 
modified as follows: (1) the COSMOS sample has been cor- 
rected for the skynoise-induced excess bias fSec. l7.2p . where 
the excess bias is tabulated in TableO (2) the modified EGS 
lens sample has been used (Sec. 17.3) to calculate the bias in 
EGS; and (3) the low resolution objects in the stripe 82 sam- 
ples (VVDS, PRIMUS-02hr and PRIMUS-23hr) have been 
removed from the source sample fSec. I3.2T3|) . Although not 
shown here, further binning of the lens objects (refiecting the 
stacking scheme, e.g. based on stellar mass or luminosity) 
would be necessary to estimate the bias for any particular 
science application. 

These panels again show that calibrating the lowest red- 
shift bin for the lens can be unreliable, as demonstrated by 
the large variance in the bias fez between the different cal- 
ibration samples. The highest redshift bin have a signifi- 
cant scatter in the bias among the subsamples which comes 
from the skynoise correction (applied to the COSMOS sam- 
ple). Otherwise, the bias correction is ~20 per cent, with 
the subsamples showing agreement to 4 per cent scatter. As 
we go to higher lens redshift bins, the efficiency and com- 
pleteness drops, but not too significantly. Purity (at 90 per 
cent) remains high in all lens redshift bins. Application of 
lens-source photo-z separation requirement does not have a 
significant effect on the lensing bias or its scatter, while com- 
pleteness (and hence efficiency) drops. This indicates that 
the minimum photo-2 separation cut as implemented here 
is unnecessary. We interpret this to mean that our weighting 
scheme fEQ. I20p sufficiently down-weights nearby lens-source 
pairs that their potentially strong calibration biases do not 
lead to significant uncertainty in the calibration bias. 

Figure [26] shows the lensing bias as a function of lens 
and source photo-z bin. The lens photo-z bins are in 0.1 
intervals, while the source photo-z calibration sample was 
split into 9 equal parts of approximately 1000 galaxies each. 
The photo-2 bin edges are given in the vertical axis of the 
figure. The highly-biased bins generally have low weighting; 
the lowest lens photo-z bin is highly biased, and carries no 
weight (hence also rendering the calibration unreliable), and 
similarly for bins where lens photo-z's are close to source 
photo-z's. This figure clearly demonstrates that the highest 
source photo-z bin carries the most weight at any lens photo- 
z bin, and also has the lowest bias (except for the lowest lens 
photo-z bin). This is expected, because as the source redshift 
is well above the lens redshift, the weight oc S"^^ increases. 



Lensing bias 



Weights 



while the the lensing strength Ej,j.;^ 
the source redshift errors. 



becomes less sensitive to 



8 SUMMARY AND CONCLUSION 

In this paper, we have sought to address technical issues 
required to use photometric redshifts for both lenses and 
sources in galaxy-galaxy lensing analyses without incurring 
significant systematic error. 

To this end, we have obtained photometric redshifts for 
SDSS DR8 5-band photometry, using the template-based 
photo-z code ZEBRA. Multiple readily available and self- 
generated photo-z catalogues were perused and compared to 
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Figure 26. Left: Lensing signal bias, as a function of photo- 
z bins in both lenses and sources. The lens bins are split into 
0.1 redshift intervals, while the source bins are split into 9 equal 
parts. Right: Lensing signal weights in each of the lens/source 
bins. The highest source photo-2 bins carry the most weight at 
any lens photo-z bin. The lowest lens photo-2 bins, as well as the 
bins where lens photo-z is near the source photo-2, are highly 
biased, but carry very little weight. 



optimize the photo-z according to the key science goals of 
the catalogue, galaxy-galaxy lensing, which led us to choose 
a method with minimal photo-z; scatter for galaxies with 
2phot > 0.1 (which dominate such analyses). The underlying 
assumption is that with a calibration sample of sufficient 
size to characterize the bias and scatter, any photo-2: biases 
can simply be calibrated out, so they should not figure into 
the choice of photo-z method (provided that they are not 
so pathological as to make the problem of calibrating them 
out with a reasonable size calibration sample intractable). 

To measure the photo-z biases and scatter, a calibration 
set was required. Tests of the calibration set, and correc- 
tions for exceptional observing conditions, were required to 
ensure that its photo-z scatters and biases accurately repre- 
sent those of the full SDSS DR8 lensing and source samples. 
The calibration sample was drawn from several effectively 
flux-limited spectroscopic samples (six for the shallower lens 
sample, five for the deeper source sample) consisting of ~9k 
galaxies overlapping the SDSS footprint. 

The measured photo-z bias and scatter allowed esti- 
mates of their contributions to bias in observed galaxy- 
galaxy lensing signals. We considered source and lens photo- 
z separately, with emphasis on the newer application to 
SDSS, lens photo-z. The use of photo-z's for source redshifts 
(in addition to lens spectroscopic redshifts) allows for grav- 
itational lensing signal calibration that is more precise than 
using an assumed redshift distribution, and allows for selec- 
tion of specific source subsamples (e.g., brighter or fainter, 
or with some lens-source separation). For a known lens red- 
shift, we estimate that the lensing tangential shear can be 
converted to surface mass density with an accuracy of 2 per 
cent. Using photo-z's for lenses will make lenses available at 
higher redshifts (out to z ~ 0.5, limited to red galaxies) and 
dimmer magnitudes (down to r ~ 21 or °"^Mr ~ —19, lim- 
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ited to blue), an increase of factors of at least ten compared 
to SDSS spectroscopic samples. Rather than lens samples 
with several tens of th ousands of galaxies per st ellar mass 
or luminosity bin as in lMandelbaum et alj (|200 6j . the sam- 
ple sizes will be of order a few xlO^, which will allow for 
very high S/N measurements (for which precise calibration 
of systematics is, therefore, a high priority). The use of lens 
photo-2 converts to surface mass density with an accuracy 
of 4 per cent when the lens redshift is 0.1 < Zphot < 0.5. 

Besides the lensing signal bias originating in distance 
relations, biases in the observable or "stacking quantities" 
(required for galaxy-galaxy lensing) arise because such quan- 
tities (absolute luminosity, rest-frame colours, stellar mass) 
also depend on redshift. These biases can be studied because 
of the availability of best-fit galaxy SED types when obtain- 
ing photo-2 from the template-based method. The SED's 
provide fc-corrections, which in turn allow for a uniform def- 
inition of absolute magnitude across a span of redshifts, as 
well as stellar mass estimates. 

Because we only have a small number of heterogeneous, 
narrow-field spectroscopic survey samples to serve as the cal- 
ibration sample, specific cautions were required. For each, 
the full depth of the photometric survey must be covered 
to accurately represent the whole sample. The wide sepa- 
ration of each of the subsamples compensates for the LSS 
in the narrow field surveys, such that the LSS can be av- 
eraged over to represent the (expected) smooth dN/dz of 
the SDSS photometric sample. We found that the differ- 
ences in survey conditions (such as seeing and sky noise) of 
the SDSS photometric sample at the locations of the spec- 
troscopic survey area can bias the magnitude completeness 
and redshift distributions for objects near the limiting mag- 
nitude (r ~ 21.8). These then affected the lensing signal 
bias, where the dim, high redshift source objects have a 
smaller bias and are highly weighted. Note that our f i nding s 
in this paper supersede those in Mandclb aum et al.l l|200a ). 
which (a) had a calibration sample that was a factor of four 
smaller, and (b) lacked corrections for the atypical observing 
conditions in SDSS at the position of the COSMOS survey. 

While the bias and scatter in photo-z or photo-z-derived 
quantities can be large, we show how calibrating them 
against the true redshifts allow us to choose a binning width 
of various lens characteristics appropriate for the uncertain- 
ties, and to estimate the correction for the bias. These biases 
can now be used to push galaxy-galaxy lensing into a new 
regime for SDSS. 



ACKNOWLEDGMENTS 

We thank Carlos Cunha, Pascal Oesch, Alex Szalay, Ist- 
van Csabai, Tamas Budavari, Jeff Newman, Feng Dong, Jim 
Gunn, Yen-Ting Lin, Michael Blanton, and Robert Feld- 
mann, for useful discussions. RN was supported in part 
by NASA LTSA grant NNG04GC90G. This work has been 
supported in part by World Class University grant R32- 
2009-000-10130-0 through the National Research Founda- 
tion, Ministry of Education, Science and Technology of Ko- 
rea. 

We thank the PRIMUS team for sharing their redshift 
catalog, and thank Alison Coil and John Moustakas for 
help with using the PRIMUS dataset. Funding for PRIMUS 



has been provided by NSF grants AST-0607701, 0908246, 
0908442, 0908354, and NASA grant 08-ADP08-0019. This 
paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 

Funding for the DEEP2 survey has been provided 
by NSF grants AST95-09298, AST-0071048, AST-0071198, 
AST-0507428, and AST-0507483 as well as NASA LTSA 
grant NNG04GC89G. Some of the data presented herein 
were obtained at the W. M. Keck Observatory, which is op- 
erated as a scientific partnership among the California In- 
stitute of Technology, the University of California and the 
National Aeronautics and Space Administration. The Obser- 
vatory was made possible by the generous financial support 
of the W. M. Keck Foundation. The DEEP2 team and Keck 
Observatory acknowledge the very significant cultural role 
and reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community and appreciate 
the opportunity to conduct observations from this mountain. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. De- 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Partic- 
ipating Institutions are the American Museum of Natu- 
ral History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve Uni- 
versity, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation 
Group, Johns Hopkins University, the Joint Institute for 
Nuclear Astrophysics, the Kavli Institute for Particle As- 
trophysics and Cosmology, the Korean Scientist Group, the 
Chinese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), 
New Mexico State University, Ohio State University, Uni- 
versity of Pittsburgh, University of Portsmouth, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 



REFERENCES 

Abazajian K. N., et al. 2009, Astrophys. J. Supp., 182, 543 
Abdalla F. B., Banerji M., Lahav O., Rashkov V., 2008, 

preprint (arXiv:0812.3831) 
Albrecht A., et al. 2006, preprint (arXiv:0609591) 
Baldauf T., Smith R. E., Seljak U., Mandelbaum R., 2010, 

Phys. Rev. D, 81, 063531 
Bartelmann M., Schneider P., 2001, Phys. Rep., 340, 291 
Bell E. F., Mcintosh D. H., Katz N., Weinberg M. D., 2003, 

Astrophys. J. Supp., 149, 289 
Benitez N., 2000, Astrophys. J., 536, 571 
Blanton M. R., Roweis S., 2007, Astron. J., 133, 734 
Brown M. J. I., et al. 2007, Astrophys. J., 654, 858 
Bruzual G, Chariot S., 2003, Mon. Not. R. Astron. Soc, 

344, 1000 
Bruzual A. G, Chariot S., 1993, Astrophys. J., 405, 538 



© 0000 RAS, MNRAS 000, 000-000 



24 Nakajima et al. 



Budavari T., Szalay A. S., Connolly A. J., Csabai I., Dick- 
inson M., 2000, Astron. J., 120, 1588 
Cacciato M., van den Bosch F. C, More S., Li R., Mo H. J., 

Yang X., 2009, Mon. Not. R. Astron. Soc, 394, 929 
Capak R, et al. 2007, Astrophys. J. Supp., 172, 99 
Coil A. L., et al. 2010, preprint (arXiv: 101 1.4307) 
Coil A. L., et al. 2004, Astrophys. J., 609, 525 
Coleman G. D., Wu C, Weedman D. W., 1980, Astro- 
phys. J. Supp., 43, 393 
CoUister A., et al. 2007, Mon. Not. R. Astron. Soc, 375, 

68 
CoUister A. A., Lahav O., 2004, Pub. Astron. Soc. Pacific, 

116, 345 
Csabai I., et al. 2003, Astron. J., 125, 580 
Cunha C. E., Lima M., Oyaizu H., Frieman J., Lin H., 2009, 

Mon. Not. R. Astron. Soc, 396, 2379 
Davis M., et al. 2003, in P. Guhathakurta ed., Proc SPIE 
Vol. 4834, Discoveries and Research Prospects from 6- and 
10-Meter-Class Telescopes IL SPIE, Bellingham, p. 161 
Davis M., Gerke B. F., Newman J. A., 2005, in Wolff S. C, 
Lauer T. R., eds, ASP Conf. Ser. Vol. 339, Observing Dark 
Energy. Astronom. Soc. Pac, San Francisco, p. 269 p. 128 
Davis M., et al. 2007, Astrophys. J. Lett., 660, LI 
Dutton A. A., et al. 2010, preprint (arXiv:1012.5859) 
Eisenstein D. J., et al. 2001, Astron. J., 122, 2267 
Eisenstein D. J., et al. 2011, preprint (arXiv:1101.1529) 
Faber S. M., et al. 2003, in lye M., Moorwood A. F. M., 
eds, Proc. SPIE Vol. 4841, Instrument Design and Per- 
formance for Optical/Infrared Ground-based Telescopes. 
SPIE, Bellingham, p. 1657 
Faber S. M., et al. 2007, Astrophys. J., 665, 265 
Feldmann R., et al. 2006, Mon. Not. R. Astron. Soc, 372, 

565 
Firth A. E., Lahav O., Somerville R. S., 2003, 

Mon. Not. R. Astron. Soc, 339, 1195 
Fu L., et al. 2008, Astron. Astrophys., 479, 9 
Fukugita M., Ichikawa T., Gunn J. E., Doi M., Shimasaku 

K., Schneider D. P., 1996, Astron. J., Ill, 1748 
Giallongo E., et al. 2005, Astrophys. J., 622, 116 
Gunn J. E., et al. 1998, Astron. J., 116, 3040 
Heymans C, et al. 2006, Mon. Not. R. Astron. Soc, 371, 

L60 
Hildebrandt H., et al. 2010, Astron. Astrophys., 523, A31 
Hirata C, Seljak U., 2003, Mon. Not. R. Astron. Soc, 343, 

459 
Hoekstra H., Hsieh B. C, Yee H. K. C, Lin H., Gladders 

M. D., 2005, Astrophys. J., 635, 73 
Hogg D. W., Baldry I. K., Blanton M. R., Eisenstein D. J., 

2002, preprint (arXiv:02 10394) 
Hogg D. W., Finkbeiner D. P., Schlegel D. J., Gunn J. E., 

2001, Astron. J., 122, 2129 
Ilbert O., et al. 2006, Astron. Astrophys., 457, 841 
Ilbert O., et al. 2009, Astrophys. J., 690, 1236 
Ivezic Z., et al. 2004, Astronomische Nachrichten, 325, 583 
Kinney A. L., et al. 1996, Astrophys. J., 467, 38 
Kleinheinrich M., et al. 2005, Astron. Astrophys., 439, 513 
Kleinheinrich M., et al. 2006, Astron. Astrophys., 455, 441 
Komatsu E., et al. 2011, Astrophys. J. Supp., 192, 18 
Le Fevre O., et al. 2005, Astron. Astrophys., 439, 845 
Leauthaud A., et al. 2011, preprint (arXiv:1104.0928) 
LeFevre O., et al. 2003, in lye M., Moorwood A. F. M., 
eds, Proc. SPIE, Vol. 4841, Instrument Design and Per- 



formance for Optical/Infrared Ground-based Telescopes. 

SPIE, Bellingham, p. 1670 
Lilly S. J., et al. 2009, Astrophys. J. Supp., 184, 218 
Lilly S. J., et al. 2007, Astrophys. J. Supp., 172, 70 
Lupton R., et al. 2001, in F. R. Harnden Jr., F. A. Primini, 

& H. E. Payne cd., ASP Conf. Ser. Vol. 238, Astronomical 

Data Analysis Software and Systems X. Astronom. Soc. 

Pac, San Francisco, p. 269 
Madau P., 1995, Astrophys. J., 441, 18 
Madgwick D. S., et al. 2003, Astrophys. J., 599, 997 
Mandelbaum R., et al. 2005, Mon. Not. R. Astron. Soc, 

361, 1287 
Mandelbaum R., Li C, Kauffmann G., White S. D. M., 

2009, Mon. Not. R. Astron. Soc, 393, 377 
Mandelbaum R., Seljak U., Baldauf T., Smith R. E., 2010, 

Mon. Not. R. Astron. Soc, 405, 2078 
Mandelbaum R., Seljak U., Hirata C. M., 2008, JCAP , 8, 

6 
Mandelbaum R., et al. 2008, Mon. Not. R. Astron. Soc, 

386, 781 
Mandelbaum R., Seljak U., Kauffmann G., Hirata C. M., 

Brinkmann J., 2006, Mon. Not. R. Astron. Soc, 368, 715 
Meiksin A., 2006, Mon. Not. R. Astron. Soc, 365, 807 
Oguri M., Takada M., 2011, Phys. Rev. D, 83, 023008 
Oyaizu H., Lima M., Cunha C. E., Lin H., Frieman J., 

Sheldon E. S., 2008, Astrophys. J., 674, 768 
Padmanabhan N., et al. 2008, Astrophys. J., 674, 1217 
Parker L. C, Hoekstra H., Hudson M. J., van Waerbeke 

L., Mellier Y., 2007, Astrophys. J., 669, 21 
Pier J. R., et al. 2003, Astron. J., 125, 1559 
Refregier A., 2003, Annu. Rev. Astron. Astrophys., 41, 645 
Reyes R., Mandelbaum R., Hirata C, Bahcall N., Seljak 

U., 2008, Mon. Not. R. Astron. Soc, 390, 1157 
Reyes R., Mandelbaum R., Seljak U., Baldauf T., Gunn 

J. E., Lombriser L., Smith R. E., 2010, Nature, 464, 256 
Richards G. T., et al. 2002, Astron. J., 123, 2945 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, Astro- 
phys. J., 500, 525 
Schrabback T., et al. 2010, Astron. Astrophys., 516, A63 
Schulz A. E., Mandelbaum R., Padmanabhan N., 2010, 

Mon. Not. R. Astron. Soc, 408, 1463 
Scoville N., et al. 2007a, Astrophys. J. Supp., 172, 1 
Scoville N., et al. 2007b, Astrophys. J. Supp., 172, 38 
Sheldon E. S., et al. 2004, Astron. J., 127, 2544 
Sheldon E. S., et al. 2009, Astrophys. J., 703, 2217 
Smith J. A., ct al. 2002, Astron. J., 123, 2121 
Stoughton C, et al. 2002, Astron. J., 123, 485 
Stoughton C, et al. 2002, Astron. J., 123, 485 
Strauss M. A., et al. 2002, Astron. J., 124, 1810 
Taniguchi Y., et al. 2007, Astrophys. J. Supp., 172, 9 
Tucker D. L., et al. 2006, Astronomische Nachrichten, 327, 

821 
Willmer C. N. A., et al. 2006, Astrophys. J., 647, 853 
Wolf C, Meisenheimer K., Rix H.-W., Borch A., Dye S., 

Kleinheinrich M., 2003, Astron. Astrophys., 401, 73 
Yoo J., Tinker J. L., Weinberg D. H., Zheng Z., Katz N., 

Dave R., 2006, Astrophys. J., 652, 26 
York D. G., et al. 2000, Astron. J., 120, 1579 
Zhang P., Liguori M., Bean R., Dodelson S., 2007, Physical 

Review Letters, 99, 141302 



© 0000 RAS, MNRAS 000, 000-000 



Lensing with photometric redshifts 25 



APPENDIX A: SELECTION PROCEDURE FOR 
CATALOGUE GENERATION 

The following are the cuts imposed on the SDSS Photo 
pipeline outputs to create the galaxy catalogues used for 
this work: 

• OBJC.TYPE = 3; 

• Require BINNED 1 in r and i, and overall; 

• Reject those with BLENDED but without NODE- 
BLEND set; 

• Reject object flags: SATURATED, SATU- 
RATED.CENTER, EDGE, LOCAL_EDGE, MAYBE.CR, 
MAYBE_EGHOST, SUBTRACTED, BRIGHT, 
TOO_LARGE, BADSKY; 

• Reject those with the following set in r and i: CR, IN- 
TER?, INTERP .CENTER; 

• Extinction-corrected rmodci < 21.8 (source catalogue) 
or 22 (full photometric catalogue; < 21 for the photometric 
lens catalogue); 

• r-band extinction Ar < 0.2; 

• Source catalogue only: Resolution factor R2 > 1/3 in 
both r and i bands (comparing the image adaptive moments 
with the PSF adaptive m oments, and requiring the object 
to be well-resolved), as in lMandelbaum et al.l (|2005l V 

Note that the cuts on the model magnitude and reso- 
lution are imposed after eliminating duplicates by choosing 
the observation with better seeing (i.e., if the object fails 
the magnitude cut in the observation with better seeing, 
and passes it in the observation with worse seeing, it will 
not be in the catalogue). 
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Figure Bl. Comparison of the bias and scatter for the "hybrid" 
and our ZEBRA photo-2 methods for the source-matched cata- 
logue. Only 50 per cent of the points are plotted for clarity. The 
limited training set in the "hybrid" method (the SDSS spectro- 
scopic samples) truncates the available photo-z's to Zp^ot < 0-6. 



APPENDIX B: 
PHOTO-Z'S 



COMPARISON WITH OTHER 



There are several other photo-a's available for the SDSS 
public catalogue. Here we illustrate the issues that made 
them sub-optimal for our galaxy-g alaxy lensing analysis . 
The publicly avai lable "hybrid" (|Abazaiian et al.l 12003 : 
ICsabai et al.l 20031), "cc2" and "dl" neural net methods 



fboth lOvaizu et al.ll200a). and t he p{z) redshift probability 
distributions ( Cunha et al.ll2009f ) are all essentially training- 
set based methods. There are a few issues related to the 
suitability of the training set for the photometric sample of 
interest; specifically, we aim for more high-z galaxies prop- 
erly classified as high- 2:, and there are simply not enough 
spectroscopic surveys available at this redshift to both train 
and test the SDSS photo-z's. 

The DR7 hybrid method uses a ~50k training set uni- 
formly selected to match the colour distribution of the pho- 
tometric sample, taken from the SDSS MAIN (r < 17.7) and 
LRG (r < 19) spectroscopic sample, including an additional 
~3k blue, high- 2: (0.25 < z < 0.4) objects to complement the 
colour distribution in the training set. This training set is, 
however, limited to bright objects, and hence to low redshifts 
(^modian ~ 0.1), which results in the generated photo-z be- 
ing similarly limited to this redshift range (see Fig. lBl|) . Our 
deeper photometric samples are demonstrably at higher red- 
shifts. Their new DR8 photo-2 sample is shown to perform 
better in the high-z regime, based on a test using half the 
spectroscopic sample to train, and the other to test (Csabai, 



private communication). We note, however, that their public 
DR8 sample is trained on the same spectroscopic sample as 
our calibration set; we require further spectroscopic redshift 
samples to independently calibrate the biases in their DR8 
sample. 

The cc2 and dl neural net methods, as well as the p(z) 
redshift probability distribution, also use the SDSS MAIN 
(435k) and LRG (>80k) spectroscopic samples as the train- 
ing set, with supplements from other spectroscopic samples 
with deeper fiux limits and broader redshift coverage, but 
smaller size (~ 7k). For the p{z) method, the training sam- 
ples are then weighted to reproduce the r-magnitude and 
color distributions of the photometric sample. The limited 
number of high-z regions contributing to photo-z estimation 
results in the LSS in those regions being imprinted onto the 
output photo-z for the full photometric sample (see Fig. lB2|) . 
The excess at z ~ 0.75 is due mostly to the non-EGS DEEP2 
spectroscopic training subsample, because they select high-z 
(z > 0.7) galaxies to target. Since the DEEP2 targeting were 
based on non-SDSS bands, the magnitude-based reweight- 
ing scheme is not complete, causing the excess structure. 
Additionally, any LSS present in the training set refiects 
directly into the p{z) (Sheldon et al. 2011, in prep.). In gen- 
eral, if the sample we use to test our calibration is also in 
the training set (which is the case here), then we expect 
excessively optimistic results. There are simply not enough 
spectroscopic samples at high-z available to split between 
testing and training. 

For the lens sample, in addition to aiming for photo-z 
accuracy, we require fe-correction for our galaxy-galaxy lens- 
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Figu re B2. Cumu lative p{z) of all SDSS galaxies, based on DR7 
p{z) jCunha et al.ll200 9V Although in the p{z) method, the cal- 
ibration samples are weighted to better represent the whole of 
the SDSS photometric sample, the LSS evident in the calibra- 
tion set appears in the p{z) photometric redshifts (cf. fig. 17 of 
ICunha et al.l 12003 for the training set distribution). The excess 
at z ~ 0.75 is due to the non-EGS DEEP2 spectroscopic training 
subsample, which intentionally targets high-z objects. (The split 
at 2 ~ 0.3 is due to an issue in the training set, and is absent in 
the DR8 version of the p{z); Sheldon et al. 2011 in prep.) 



ing analysis, to stack lenses according to fc-corrected galaxy 
observables such as luminosity and colour. Such information 
is not available with the dl and cc2 methods, despite hav- 
ing the least bias and scatter for the lens sample. Although 
training-set bas ed method can, in principle, perform spectra l 
classifications (|Firth et all l2003l : ICoUister fc Lahli^ |2004 ). 
the estimated SED it is not available with these photo-z's, 
and hence fc-correction is not available. The hybrid method 
includes an SED estimator which enables a fc-correction, but 
their Zphot scatter at low photo-z is higher than with the ZE- 
BRA method. 

We briefly compare our photo-z to the MegaZ- 
LRG sample by ICollister et al.l (|2007l ) with a scatter of 
<^Az/{i+z) ~ 0.045. We impose the MegaZ-LRG color and 
magnitude cuts to our calibration sample, and find that ZE- 
BRA classifies the calibration sample into Ell and Sbc types 
with 259 and 102 objects, with scatter of crAz/(i-|-z) = 0.04 
and 0.07, respectively. For the combined sample, the scat- 
ter is 0.05. While the scatter in our calibration sample is 
not quite as small as that from MegaZ-LRG, our cuts do 
not include the fiber magnitude cut or the magnitude-based 
star-galaxy separation criteria cuts imposed on the MegaZ- 
LRG sample, which presumably would decrease the number 
of outliers we see in our photo-z error. Additionally, we note 
that our need for fc- correct ions to define luminosities pre- 
vents us from using photo-z estimators witho ut any template 
inform ation, such as the neural net method of lCollister et al.l 
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